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FOREWORD 

This  report,  the  fourth  in  the  Boulder  Atmospheric  Observatory 
(BAO)  series,  brings  together  in  one  volume  five  studies  of  the  noc- 
turnal stable  layer.   The  work  was  conducted  in  the  Wave  Propagation 
Laboratory  during  1981  and  1982.   Impetus  for  the  first  four  of  these 
studies  came  from  visiting  scientists  who  found  in  data  collected  at 
BAO  a  unique  source  of  information  on  mean  and  turbulent  character- 
istics of  stable  layers.   The  structure  of  turbulence  in  nocturnal 
stable  layers,  unlike  that  of  daytime  convective  boundary  layers, 
cannot  be  described  in  simple  universal  terms.   Stable  layers  are 
almost  always  in  a  state  of  evolution.   At  BAO,  the  problem  is  even 
more  challenging  because  of  the  rolling  terrain,  the  drainage  flow 
from  the  Rocky  Mountains,  and  the  frequent  presence  of  gravity  waves 
and  elevated  inversion  layers.   To  the  scientist  eager  to  unravel 
these  complexities,  the  BAO  records  offer  conditions  of  unending 
variety.   Each  paper  in  this  report  explores  a  different  aspect  of 
the  nocturnal  stable  layer.   It  is  our  hope  that  this  collected 
work  will  stimulate  comments  from  other  workers  in  the  field  and 
lead  to  new  insights  into  the  structure  of  the  stable  atmosphere. 


J.C.  Kaimal,  Editor 

Atmospheric  Studies  Program  Area 

Wave  Propagation  Laboratory 
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1.   OBSERVATIONS  OF  TURBULENCE  STRUCTURE  IN  STABLE  LAYERS 
AT  THE  BOULDER  ATMOSPHERIC  OBSERVATORY 

J.C.R.  Hunt 
Department  of  Applied  Mathematics 

and  Theoretical  Physics 

University  of  Cambridge 
Cambridge  CB39EW,  United  Kingdom 

J.C.  Kaimal,  J.E.  Gaynor,  and  A.  Korrell 

NOAA/ERL/Wave  Propagation  Laboratory 

Boulder,  Colorado   80303 


ABSTRACT.   Three  days  of  measurements  in  stable  conditions  (Monin- 
Obukhov  length  L  between  15  and  300  m)  at  the  Boulder  Atmospheric 
Observatory  (BAO)  tower  are  presented.   Winds  came  off  the  Rocky 
Mountains  on  2  days  and  off  the  plains  on  the  other  day.   Vertical 
profiles  of  the  mean  horizontal  velocity  U,  and  the  Brunt-Vaisala 
frequency  N,  and  of  the  variances,  fluxes,  spectra,  cospectra, 
quadspectra,  and  correlations  are  examined  in  this  study.   New 
theoretical  arguments  are  developed  for  relating  variances  to 
scales  of  temperature  6  and  vertical  velocity  w  fluctuations  and  to 
scales  of  heat  and  momentum  fluxes  in  the  stable  boundary  layer. 
Some  tentative  conclusions  drawn  from  this  study  are  the  following: 

(1)  In  terrain,  even  gently  rolling_as  it  is  at  BAO 
(slopes  <7%),  the  profiles  of  dU/dz  and  d6/dz  are  significantly 
different  from  those  over  level  ground.   For  example,  in  stable 
conditions,  z(dU/dz)  and  z(d6/dz)  often  do  not  increase  with 
height  z. 

(2)  The  presence  of  coherent  wave  motions  cannot  be  detected 
from  U  and  6  profiles.   The  most  sensitive  tests  are  w0  cospectra 
and  quadspectra. 

(3)  Wave  motions  may  be  apparent  only  in  the  9  spectra  (mod- 
erate wave  conditions)  or  may  be  apparent  in  both  w  and  8  spectra 
(strong  wave  conditions).   Typical  inertial  subrange  spectra  are 
found  in  all  cases. 

(4)  Whether  or  not  waves  are  present,  a  theoretical  prediction 
for  local  parameterization  of  temperature  fluctuation  Qfi  and  heat 
flux  w0  based  on  O    ,  N,  and  d6/dz  appears  to  order  the  data  fairly 
well  where  strong  uniform  density  gradients  exist.   In  particular 

we  find  the  dimensionless  temperature  fluctuation  parameter  Co  ~ 
0.8  ±  0.25  and  thermal  diffusivity  parameter  F  -  0.2  ±  0.1. 


(5)  In  the  boundary  layer,  where  dU/dz  >  N,  the  integral 
scale  of  the  w  component  L    is  found  to  agree  (within  about  50%) 

X 

with  the  theoretical  formula, 

[L(w)-j-l  a  Q7    (du/dz)/a  +  z-1. 
X  w 

(6)  The  large  values  of  the  cospectra  of  w  and  6  at  low 
frequencies  show  that  wavelike  motions  can  transfer  significant 
quantities  of  heat  (or  pollutant).   This  is  probably  due  to 
small-scale  mixing  on  a  time  scale  of  about  5  N   accompanying 
the  wave- induced,  large-scale  motions  across  the  temperature 
gradient. 


1.1   INTRODUCTION 

The  general  aim  of  this  paper  is  to  present  some  interesting  measure- 
ments, taken  on  the  Boulder  Atmospheric  Observatory  (BAO)  300-m  tower,  of  the 
mean  and  fluctuating  velocity  and  temperature  in  stably  stratified  atmospheric 
air  flows  within  300  m  of  the  ground.   This  depth  may  or  may  not  exceed  the 
depth  of  the  stable  boundary  layer.   From  the  analysis  of  these  measurements, 
a  number  of  general  features  of  these  flows  can  be  discerned;  we  particularly 
emphasize  those  features  that  are  of  practical  importance  for  estimating  flow 
and  diffusion  in  stable  conditions,  both  over  level  ground  and  in  hilly  ter- 
rain, and  those  features  of  the  flow  that  are  affected  by  the  terrain  at  BAO. 
The  BAO  tower  is  situated  in  rolling  terrain  30  km  from  the  foothills  of  the 
Rocky  Mountains.   It  is  of  interest  to  know  what  features  of  the  flow  are 
similar  to  those  measured  over  the  flat  terrain  at  Kansas,  Minnesota,  and  St. 
Louis,  Mo.   The  results  presented  here  are  not  chosen  with  the  view  to  study 
the  development  of  the  boundary  layer;  rather  we  concentrate  on  the  structure 
of  these  flows  averaged  over  periods  on  the  order  of  20  or  60  min. 

The  measurements  bring  out  some  interesting  features  of  the  large-scale 
turbulence  in  and  above  the  stable  boundary  layer  (SBL).   In  the  upper  part 
of  the  SBL,  in  all  kinds  of  terrain,  there  appears  to  be  a  significant  amount 
of  energy  in  internal  wave  motion  mixed  in  with  the  turbulence.   By  use  of 
various  statistical  analyses  we  attempt  to  identify  these  various  kinds  of 
fluctuations,  to  estimate  the  relative  energies  in  these  fluctuations,  and 
to  find  out  their  relative  contributions  to  the  heat  fluxes  and  to  the  de- 
struction of  turbulent  energy  by  buoyancy  flux  and  by  viscous  dissipation. 


More  than  in  convectively  unstable  or  neutral  flows,  the  dynamics  of  a 

stably  stratified  turbulent  flow  depend  on  the  rate  at  which  fluid  elements 

mix  with  each  other  and  change  their  density.   If  there  is  no  such  mixing, 

the  vertical  displacements  of  fluid  elements  are  constrained  by  the  amount  of 

energy  in 'the  turbulence  to  lie  within  a  distance  on  the  order  of  a  /N  of 

w 

their  equilibrium  static  height,  where  a  is  the  standard  deviation  of  ver- 
tical velocity  and  N  is  the  Brunt-Vaisala  frequency  (Pearson  et  al. ,  1983; 
Pearson  and  Britter,  1980).   By  measuring  the  vertical  velocity  distribution 
and  spectra,  the  temperature  fluctuations,  and  the  heat  flux,  and  by  making 
use  of  the  Lagrangian  model  of  Csanady  (1964)  and  Pearson  et  al. ,  we  can 
learn  a  good  deal  about  the  particle  displacements  and  mixing.   We  are  inter- 
ested in  knowing  to  what  extent  internal  wave  motion  affects  these  relations. 
From  these  measurements,  predictions  (by  means  of  the  Lagrangian  model)  for 
vertical  diffusion  of  plumes  can  be  made,  which  we  hope  to  test  later  in  the 
forthcoming  BAO  diffusion  experiments. 

Another  aspect  of  the  turbulence  structure  we  measure  and  analyze  here 
is  the  effect  of  stable  stratification  on  the  length  scales  of  velocity  and 
temperature.   Using  an  analysis  (unpublished  but  summarized  in  Appendix  B) 
based  on  the  turbulence  kinetic  equation,  we  examine  the  main  effects  of 
stable  stratification  on  the  turbulence  in  the  lower  part  of  the  boundary. 
Since  the  length  scales  affect  the  dissipation  of  turbulent  energy  and  the 
destruction  of  temperature  variance,  there  are  some  important  dynamical  and 
modeling  consequences  from  our  findings.   Inter  alia,  we  show  that,  in  the 
lower  part  of  the  stable  boundary  layer,  the  main  effect  of  an  increase  in 
stable  stratification  on  the  turbulent  energy  is  not  so  much  the  loss  of 
kinetic  energy  by  the  buoyancy  flux  (the  usual  explanation)  as  the  larger  (by 
a  factor  of  4)  increase  of  viscous  dissipation  caused  by  the  reduction  in 
length  scale  associated  with  the  increased  mean  velocity  gradient. 


1.2   TOPOGRAPHIC  AND  SYNOPTIC  CONDITIONS  OF  THE  EXPERIMENTS 


The  measurements  were  made  on  the  BAO  tower  about  30  km  east  of  the 


foothills,  and  about  60  km  from  the  Continental  Divide  of  the  Rocky  Mountains, 
where  the  mountains  are  at  an  effective  height  H  (-1.5  x  10  m)  above  the 
ground  level  at  BAO.   In  stable  conditions,  unlike  unstable  conditions,  the 
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Figure  1.1. — Topography  of  the  BAO  site  and  surroundings:  (a)  E-W  sec- 
tion through  the  BAO  site;  (b)  map  showing  the  location  of  BAO  with 
respect  to  the  Continental  Divide. 


effects  of  large  hills  can  affect  the  flow  many  kilometers  downwind  or  upwind; 
also,  the  flow  over  the  local  topography  (where  slopes  are  at  most  of  order 
7%)  has  to  be  considered  more  carefully  (Kaimal  et  al. ,  1982).   Details  of 
the  topography  are  given  in  Figs.  1.1  and  1.2. 

Surface  weather  maps  and  rawinsonde  plots  bracketing  the  observation 
periods  on  18  April  1978,  22  April  1978,  and  15  April  1980,  and  upper-air 
soundings  at  Denver  (30  km  from  BAO)  are  given  in  Appendix  A. 


On  the  second  and  third  days  analyzed,  22  April  1978  and  15  April  1980, 

the  general  air  flow  over  300  km  around  the  site,  and  the  local  air  flow, 

were  from  the  west  to  northwest  and  were  therefore  coming  from  the  mountains. 

The  tower  is  at  a  distance  of  about  30  to  40  times  H  from  the  mountains. 
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Figure  1.2. — A  conventional  contour  map  of  the  immediate  BAO 
(elevation  in  feet). 


terrain 


Laboratory  experiments  (e.g.,  Hunt  and  Snyder,  1980),  calculations  (e.g., 
Klemp  and  Lilly,  1975),  and  aircraft  measurements  show  that  the  effects  of 
the  Rockies  on  the  wind  extend  at  least  as  far  as  this  in  stable  conditions. 
So  on  these  days  we  expect  the  air  flow  at  BAO  to  have  characteristics  quite 
different  from  that  over  a  plain,  whereas  in  the  unstable  conditions  measured 
by  Kaimal  et  al.  (1982)  such  winds  behaved  much  as  they  did  over  the  plains 
of  Kansas  or  Minnesota.   [It  would  be  interesting  to  analyze  the  upper-air 
soundings  given  in  Appendix  A  to  see  if  and  what  kinds  of  internal  gravity 
waves  are  expected;  perhaps  calculations  such  as  those  performed  by  Einaudi 
and  Finnigan  (1981)  might  be  possible.] 


1.3  MEAN  WIND  SPEED  AND  STRATIFICATION  DATA 

Figure  1.3  shows  the  vertical  profiles  of  the  mean  horizontal  velocity  U 

and  direction  A,  and  the  mean  potential  temperature  gradient  d6/dz  calculated 

_        1/2 

in  the  form  of  the  Brunt-Vaisala  frequency  N  =  [g(d0/dz)/T  ]    ,  where  T  is 

the  ground- level  temperature  and  g  is  the  acceleration  due  to  gravity.   On 
the  first  day  of  measurements,  18  April  1978,  the  wind  came  from  between 
north  and  northeast  both  locally  and  generally.   The  only  reason  in  this  case 
to  expect  results  different  from  those  over  a  plain  is  because  of  the  local 
terrain  around  the  tower,  which,  as  we  shall  show,  is  sufficiently  sloping  to 
have  a  measurable  effect  in  stable  conditions.   On  22  April  1978  and  15  April 
1980  the  wind  came  off  the  mountains.   (Numerical  values  of  U  and  other  mean 
and  fluctuating  velocity  and  temperature  data  are  presented  in  Appendix  A.) 

A  striking  feature  of  the  profiles  of  U  is  how  only  in  the  third  case, 
where  there  was  a  collapsing  daytime  mixed  layer,  is  the  mean  velocity  shear 
dU/dz  effectively  confined  to  about  200  m.   In  the  second  case,  in  the  middle 
of  the  night  with  a  strong  wind,  there  is  appreciable  shear  all  the  way  up 
to  300  m.   (The  value  of  L,  the  Monin-Obukhov  length,  is  about  the  same  at 
1700  MST  on  the  first  day  as  between  0006  and  0436  MST  on  the  second  day).   In 
the  first  case  considered  here,  18  April  1978,  the  shear  is  significant  at  the 
top  of  the  tower.   Only  the  profiles  on  15  April  1980  are  similar  to  the 
Minnesota  profiles  measured  by  Caughey  et  al.  (1979),  where  dU/dz  was  usually 
negligible  above  about  100  m  in  strongly  stable  conditions. 

The  upper-air  soundings  at  1700  on  18  April  and  at  0500  on  19  April 
1978  were  too  different  to  infer  the  conditions  at  2100  to  2300.   On  22  April 
1978  there  was  also  little  change  in  the  wind's  direction  or  strength,  but  the 
value  of  N  varied  with  height  z:  N  -  0.007  rad/s  for  z  <  0.8  km;  0.001  rad/s 
for  0.8  km  <  z  <  1.5  km;  and  0.037  rad/s  for  1.5  km  <  z  <  3  km.   On  15  April 
1980  the  soundings  show  that  the  direction  and  strength  of  the  wind  and  the 
value  of  N  did  not  change  significantly  up  to  about  4  km  above  the  ground. 

To  what  extent  are  the  tower  profiles  affected  by  the  local  terrain  when 
the  winds  are  from  the  west,  given  a  local  upslope  of  about  2%  to  the  west  of 
the  tower?   The  effect  of  a  slope  on  U(z)  depends  on  the  shear  in  the  approach 


t>  O  <J 


o 
co 

l-  I-  H 
W  W  CO 

O) 

552 

o  o  o 

owo 

< 

S00O1 

1 

>  o  < 

1 

o 

o 

o 

o 

en 

C\J 

(Ul)Z 


(UJ)Z 


>  o< 


in 
o 
o 

o 

CO 

■6 

X) 
0) 

ra 

en 

w 

co 

Z 

to 

o 

H 

•H 

o 

o 

O 

CO 

J-t 

Cu 

o 

CO 

en 
0 

•■> 

cvj 

!Z3 

< 

o 

CD 
CM 

o      • 

c  o 

<U    00 
CT   rH 

CO 

o 

M     rH 

yancy  f 
15  Apri 

CO 

o 

m 

E 

3  /-s 

,£)    X) 

Z> 

)  ,  and 
1978;  ( 

o 

9 

CO 

o 

d 

direction  (A 
(c)  22  April 

CM 

^ 

XI 

c_> 

^ 

C    XI 

CD 

•h  a 

— 

S     cfl 

XI 

00 

0) 

r~» 

a) 

CT. 

a,  r-\ 

CO 

H 

X) 

•H 

fi 

M 

•H 

fX 

5 

<U 

H 

00 

CO 

rH 

4J 

C 

/ — V 

o 

CO 

N 

■> — ' 

•H 

5-1 

•  • 

O 

CO 

IS 

a) 

1 

&D 

1 

CO 

• 

5-i 

m 

CI) 

* 

> 

H 

cd 

0) 

a 

M 

■H 

3 

0 

ox) 

1 

•H 

O 

Pn 

CN 

(LU)Z 


flow.   In  stable  conditions  the  value  of  the  velocity  shear  in  the  approach 
flow  is  far  greater  than  its  value  in  neutral  conditions.   Calculations  and 
wind  tunnel  studies  indicate  that  U(300) /U(10)  is  typically  50%  greater  in 
stable  conditions  than  in  neutral  ones  (Hunt,  1981).   This  can  increase  the 
local  normalized  speedup  of  the  mean  wind  speed  AU/U  by  a  factor  of  2  or  more. 
One  might  expect  AU/U (10  m)  to  be  about  0.02  [U(300)/U(10) ]  -   0.2.   The  dy- 
namical effects  of  the  density  gradient  on  the  flow  over  this  slope  are  not 
large  since  U(300)/ [N(300)  x  L  ]  >  1  (Hunt,  1981).   L  is  the  half-length  of 
the  equivalent  hill  (-800  m  in  this  case).   This  increase  in  the  speedup 
reduces  the  shear  over  a  hill  or  slope  and  perhaps  explains  why  on  18  and 
22  April  1978  z(dU/dz)  increases  more  slowly  with  z  than  is  observed  on 
level  ground.   The  plot  of  kz(dU/dz)/u^  as  a  function  of  z/L  in  Fig.  1.4a 
shows  an  increase  with  z,  as  is  found  over  level  ground,  but  with  a 
smaller  coefficient:  a  =  3.2  rather  than  5.   Here  k  is  the  von  Karman's 
constant  and  u^  is  the  surface  friction  velocity.   Our  measurements  of 
z(dG/dz)  do  not  increase  with  height,  unlike  the  measurements  over  level 
ground  (Kaimal,  1973). 

An  important  parameter  for  the  calculation  of  air  flow  over  hills  is 
the  Scorer  (1949)  parameter  £~2  =  [(N/U)2  -  (d2U/dz2)/U] .   In  many  calcu- 
lations it  is  assumed  that  in  stable  conditions  I     is  largely  determined 
by  (N/U)  for  z  >  h,  h  being  the  thickness  of  the  stable  surface  layer  (e.g., 
Smith,  1980). 

2        2    2 
The  calculations  of  (N/U)   and  (-d  U/dz  )/U  are  plotted  in  Fig.  1.4b 

and  show  that,  for  the  3  days  considered  here,  the  mean  velocity  curvature 
largely  determines  I     for  z  <  300  m.   Consequently  in  rolling  or  hilly  ter- 
rain, when  calculations  are  to  be  made  of  air  flow  over  a  specific  hill  or 
over  a  limited  region  of  hills,  the  effects  of  the  incident  shear  profile  may 
well  be  as  important  as  the  dynamical  effects  of  the  stable  stratification. 
The  curvature  may  extend  farther  up  into  the  boundary  layer  in  stable  con- 
ditions than  in  neutral  conditions. 
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Figure  1.4. — Analysis  of  mean  velocity  and  density  profiles:  (a)  com- 
parison of  the  dimensionless  velocity  gradient  with  the  Monin-Obukhov 
similarity  law  for  18  April  1978  (k  =  0.4  is  assumed);  (b)  comparison 
of  the  profile  curvature  term  (-d2U/dz2)/U  term  with  (N/U)2. 


1.4   ANALYSIS  OF  VELOCITY  AND  TEMPERATURE  FLUCTUATIONS  AND  HEAT  FLUX 
IN  TERMS  OF  WAVES  AND/OR  TURBULENCE 

A  cursory  look  at  Figs.  1.5-1.8  shows  that  the  fluctuations  in  vertical 
velocity  and  temperature  exhibit  waves  as  well  as  turbulence,  waves  being 
defined  in  the  sense  of  regular  oscillations.   (Since  these  were  tower  mea- 
surements, we  could  not  make  a  more  decisive  test  of  waves  by  distinguishing 
whether  these  oscillations  were  or  were  not  propagating  relative  to  the  mean 
flow.)   On  the  days  chosen  here  for  analysis,  varying  amounts  of  wave  motion 
relative  to  the  turbulence  are  observed;  the  wave  conditions  are  weak,  moder- 
ate, or  strong  on  these  occasions,  and  are  classified  accordingly. 

In  Fig.  1.5a,  vertical  velocity  w  and  temperature  9  fluctuations  at  10 
and  150  m  are  plotted  as  functions  of  time  for  22  April  1978.   Note  the  strong 
strong  wavelike  behavior  in  w  at  150  m.   By  contrast,  in  Fig.  1.5b  the  traces 
for  15  April  1980  show  weak  wavelike  behavior.   The  waves  are  certainly  or- 
ganized, but  their  amplitude  is  weak. 

In  Figs.  1.6a,b  the  spectra  of  w  and  0 — nS  (n)  and  nS„(n) — and  the  co- 

W  u 

spectra  and  quadspectra  of  w9 — nC  ft(n)  and  nO  fl(n) — are  plotted  for  the 
turbulence  on  18  April  1978.   Note  that  there  appear  to  be  concentrations 
of  energy  at  discrete  frequencies  in  all  these  spectra,  a  fact  that  can  be 
tested  by  looking  for  discrete  frequencies  in  autocorrelation  functions. 
There  is  a  much  greater  contribution  by  these  wavelike  fluctuations  to  the 
temperature  than  to  the  velocity  spectra.   The  object  of  plotting  both  the 
cospectra  and  the  quadspectra  is  to  show  how  much  heat  flux  is  actually 
transferred  by  the  low- frequency,  wavelike  motions.   For  a  nearly  pure  sine 
wave  with  no  mixing  in  a  temperature  gradient,  the  cospectra  would  be  very 
small  compared  with  the  quadspectra.   In  Figs.  1.6a,b  we  observe  that,  com- 
pared with  the  velocity  spectra,  a  significant  portion  of  the  heat  flux  is 
transferred  by  the  low- frequency ,  wavelike  motions,  because  C  o(n)  is  of  the 


same  order  as  0  fi(n).   It  is  also  observed  that,  at  the  frequencies  where 

— 2 
the  6   "energy"  is  large*  and  Sfl(n)  is  near  its  maximum,  0  „(n)  has  the  foi 

of  a  double  "spike"  as  it  rapidly  swings  from  one  sign  to  the  other  over  a 


*  The  drift  in  6  produces  about  10%  of  the  variance, 
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narrow  frequency  range.   This  "spiky"  form  is  indicative  of  large  phase  shifts 
between  w  and  6  and  is  a  nonlinear  phenomenon. 

The  w  spectrum  in  Fig.  1.6c,  for  15  April  1980,  indicates  that  on  this 
occasion,  when  L,  N,  and  U  have  about  the  same  values  as  on  18  April  1978, 
there  is  no  discernible  wave  motion.   But  the  6  spectrum  in  Fig.  1.6c  and 
the  w6  cospectrum  and  quadspectrum  in  Fig.  1.6d  indicate  some  weak  concen- 
tration of  energy  into  bands  of  frequencies.   Note  that  0  fi(n)  does  not  have 
a  spiky  form  at  frequencies  where  C  n(n)  i-s  a  maximum,  another  indication 
of  the  weakness  of  wave  motion.   For  frequencies  where  the  wave  motions  are 

weak,  it  is  possible  to  define  a  length  scale  L  =  U/2iTn  even  for  the  tem- 

x        m 

perature  fluctuations.   Here  n  is  the  value  of  frequency  n  where  S(n)  or 

m 

C(n)  is  a  maximum.  The  turbulent  flow  with  very  weak  waves  is  comparable  with 

the  Kansas  air  flow  analyzed  by  Kaimal  (1973).   The  graphs  are  not  shown 

here.   However,  it  was  found  that  the  ratio  L   /L    is  in  the  range  6  to  8 

xx 

at  22  and  50  m,  but  this  ratio  is  only  about  3  at  150  m  (as  may  be  inferred 
from  Fig.  1.6c);  l/W  Vl(w)  is  about  2  at  150  m. 


Figure  1.7  shows  typical  spectra,  cospectra,  and  autocorrelations 

R  (t)  for  strong  wave  conditions.   Figure  1.7a  shows  how  the  strength  of 
w 

the  wave  motion  centered  at  a  few  frequencies  increases  with  height.   (These 
frequencies  can  be  picked  out  from  the  autocorrelations  in  Figs.  1.7c,d.) 
At  10  m  the  w  spectrum  has  much  the  same  form  as  the  moderate  and  weak  wave 
spectra  of  Figs.  1.6a,c.   At  greater  heights,  note  the  concentration  of  the 
heat  flux  at  a  few  frequencies  (Fig.  1.7b)  and  the  spiky  shape  of  the  Q  fl(n) 
quadspectrum.   The  contribution  of  these  wavelike  motions  to  w6  is  even 
greater  than  in  the  moderate  wave  conditions  shown  in  Fig.  1.6b. 

Figures  1.8a,b  show  how  the  energy  in  the  cospectrum  and  quadspectrum  of 
the  vertical  velocity  measured  at  100  and  150  m  is  also  concentrated  around  a 

few  frequencies,  the  first  and  second  most  energetic  waves  (denoted  by  n  and 

-3  -2  * 

n_)  being  at  about  7  x  10   and  5  x  10   Hz  respectively.    In  Fig.  1.8c,  we 

explore  the  variation  of  the  energy  in  these  two  modes  with  height.   We  ob- 
serve a  noticeable  maximum  in  nS(n=n..)  at  z  =  200  m,  but  not  in  nS(n=n?). 

*  The  first  of  these  frequencies  is  close  to  the  value  of  N/2tt  in  the  deep 
stratified  layer  between  1.5  and  3  km. 

15 


The  latter,  monotonic  increase  is  somewhat  similar  to  the  energy  in  the  wave 
studied  by  Einaudi  and  Finnigan  (1981).   The  former's  distribution  is  distinctly 
dissimilar.   The  correlation  of  the  vertical  turbulence  and  temperature,  and 
also  the  correlation  spectra  of  vertical  velocity  at  the  first  two  most 
energetic  frequencies,  are  plotted  in  Fig.  1.8d.   Note  how  the  correlation 
coefficient,  R     (subscript  4  denotes  100-m  level),  falls  off  with  I z— 100  i , 

WW.  ll' 

z  4 
the  height  above  or  below  100  m,  i.e.,  at  a  rate  similar  to  the  falloff  of 

C    (n=n„).   But  the  correlation  in  C    (n=nn )  is  much  higher  and  is  com- 
w  w.     2  w  w.     1  ° 

z  4  z  4 

parable  with  that  of  L  „  ,  showing  again  how  the  temperature  fluctuations 

z  4 

are  driven  by  wave  motions  to  a  greater  extent  than  the  vertical  velocity 

fluctuations  are. 


Estimating  the  relative  amounts  of  energy  in  turbulence  and  waves  cannot 
be  unambiguous  since  waves  may  well  have  random  phases  and  be  rotational,  two 
features  that  are  also  important  in  low-frequency  turbulence.  If  waves  are 
motions  with  energy  centered  over  a  narrow  frequency  range,  they  can  be  identi- 
fied if  the  cospectra  and  quadspectra  are  spiky.   To  estimate  the  energy  of 
turbulence  in  the  presence  of  waves,  we  assume  that  the  contribution  of  tur- 
bulence to  the  spectrum  is  of  the  same  form  as  that  in  the  absence  of  waves 
(Figs.  1.9a,b).   This  was  the  assumption  of  Caughey  et  al.  (1979)  in  their 
analysis  of  the  Minnesota  data.   Thence  we  can  say  that  on  the  basis  of 

Figs.  1.6a,c  on  18  April  1978  and  15  April  1980  for  z  <   300  m,  and  of  Fig.  1.7a 

2 
on  22  April  1978  for  z  <  10  m,  the  contribution  of  waves  to  the  variance  0 

~  2 

was  weak.   But  waves  provided  about  50%  of  O     at  150  m  on  the  strong-wave 

occasion  of  22  April  1978. 

A  further  distinction  between  different  kinds  of  wave-turbulence  com- 
binations appears  on  examining  the  temperature  and  heat  flux  spectra.   As 
is  shown  schematically  in  Figs.  1.9a-d  on  different  stable  days,  w  spectra 
may  be  largely  similar  to  those  in  turbulent  shear  flows,  but  0  spectra  and  w9 
cospectra  may  be  quite  different:  on  one  occasion  dominated  by  waves  and  on 
the  other  not.   In  the  first  of  these  two  situations,  we  shall  refer  to  the 
waves  as  "strong"  (22  April  1978)  or  "moderate"  (18  April  1978);  on  the 
second,  we  shall  refer  to  them  as  "weak"  (15  April  1980). 
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Figure  1.9. — Schematic  form  of  spectra  in  stable  conditions:  (a)  vertical 
velocity  in  the  presence  of  weak  waves  (15  April  1980);  (b)  vertical 
velocity  with  strong  waves  (22  April  1978);  (c)  heat  flux  spectra;  (d) 
temperature  spectra. 


To  understand  the  dynamical  processes  in  a  stratified  turbulent  flow,  one 
must  estimate  the  rate  of  energy  dissipation  £  per  unit  mass,  and  relate  it 
to  the  rate  of  buoyant  production  per  unit  mass,  B  =  g(w8)/T  ,  and  to  the 

ratio  of  the  macroscopic  velocity  and  integral  scales,  0.A  a  /L  '  (see 

Appendix  C).   It  is  also  of  interest  to  estimate  which  part  of  the  energy 

spectrum  is  contributing  most  to  B.   Only  if  the  major  contribution  to  B 

occurs  at  frequencies  less  than  n  can  £  be  estimated  from  the  inertial  sub- 
in 

range  of  the  spectrum.   (Otherwise  the  turbulent  energy  in  the  subrange  would 
be  transforming  itself  into  potential  energy.) 


Our  analysis  of  the  measurements  of  C  qC11)  shows  that,  even  in  the 
weak-wave  case  where  L 


(w6) 


(w) 
2L    ,  the  turbulence  with  scales  larger  than 


(w) 


produces  almost  all  the  buoyancy  flux.   Consequently  £  should  be  given 


by  the  value  of  S  (n)  in  the  inertial  subrange  according  to  the  result 

S  (n)  -  (4/3)  0.5  £2/3(2TTn/U)'5/3  (Kaimal  et  al.,  1972). 
w 
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Using  our  measurements  of  18  April  1978  at  2120  MST  (moderate  waves) , 

2   3  3  (v)  2      3 

we  estimate  £  -  0.0046  m  /s  at  50  m  and  0.4  a   /Lv  '-   0.005  m  /s  .   For  this 

w  x 

case  we  find  | B/e |  -  1/5.   On  the  strong-wave  occasion  of  22  April  1978  at 

150  m,  e  -  0.0037  m  /s  and  | B/e |  -  1/2.   It  is  interesting  to  compare 

3  -1 
(e/0.4  a  )   with  two  length  scales:  the  scale  of  the  peak  of  the  nS  (n) 

spectrum,  2frn  /U,  and  the  scale  0  /N.   In  this  case,  inspection  of  the 

spectra  in  Fig.  1.7  shows  that  the  turbulent  component  of  the  w  variance  is 

3     -1  ~ 

about  half  the  total  variance.   Thence  the  scale  [e/(0.4  a  .  s)]    =  59  m, 

-1  -1  wvt; 

whereas  (27Tn  /U)  "  -  40  m  and  (N/a  )   -  100  m.   They  are  all  of  the  same 
m  w 

order  (as  they  are  in  non-wavy  conditions). 


Since  in  strongly  stable  conditions  the  integral  scale  of  the  vertical 
component  of  turbulence  is  small  compared  with  the  distance  above  the  surface 
(for  z  >  25  m) ,  it  is  plausible  that  the  turbulence  scale  is  largely  deter- 
mined by  local  conditions,  namely  the  values  of  dU/dz,  O  ,  and  N,  which 

w 

determine  the  two  natural  inverse  length  scales  (dU/dz)/a  and  N/a   (see 

W(w)     W 
Fig.  1.9).   (Theoretical  arguments  for  the  dependence  of  L    on  z,  a  ,  dU/dz, 

and  N  are  summarized  in  Appendix  D.)   An  integral  scale  is  most  easily  de- 

(w) 

rived  by  using  the  spectrum  of  w  to  infer  that  L    -  U/(2fTn  ),  where  n  is 
&      r  x  m         m 

the  value  of  n  where  nS  (n)  is  a  maximum.   (This  definition  is  equal  to  the 

w 

true  integral  scale  if  the  autocorrelation  is  an  exponential  function.)   Evi- 
dently this  definition  is  not  appropriate  for  turbulence  with  strong  waves, 
where  the  spectrum  (as  in  Fig.  1.7a)  does  not  have  a  maximum. 


The  values  of  n  ,  even  in  the  spectra  chosen  here,  are  uncertain  to 
m'  r  ,  . 

(w) 
within  a  factor  of  about  +25%.   First,  Figs.  l.lOa-c  show  how  L    increases 

—  x 

to  a  large  but  finite  value  where  the  shear  vanishes  (dU/dz  ■*  0).   This  finite 

value  approximates  to  a  /N.   Second,  the  results  show  that  in  the  lower  part 

W  (w)  -1 

of  the  shear  layer  the  vertical  profile  of  [L   ]   more  closely  approximates 

the  profile  of  (dU/dz)/a  than  the  profile  of  N/a.   To  within  50%,  the  BAO 

w 

data  agrees  with  the  formula  for  level  ground  [(D.14)  of  Appendix  D]  that 

[L(w)]_1  -  0.7  (dU/dz)/a  +  1/z  (when  dU/dz  >  N) . 
x  w 

(w) 
Note  in  Fig.  1.10c  how  the  values  of  L    measured  in  the  Kansas  ex- 

x 

periment  at  the  same  Monin-Obukhov  length  L  are  about  half  the  values  of 
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L    in  these  experiments,  where  wave  motions  (even  though  weak)  may  increase 
L   .   Despite  the  scatter  in  our  results,  this  difference  is  probably  real. 

The  spectra  of  6  and  cospectra  of  w0  in  Figs.  1.6  and  1.7  show  how  the 
temperature  fluctuations  and  the  heat  flux  are  controlled  by  motions  with 
larger  scales  than  that  of  the  vertical  velocity.   The  results  of  Fig.  1.10 
indicate  that  such  scales  are  likely  to  be  affected  by  the  buoyancy  length 

scale  o  /N.   In  the  moderate-  and  strong-wave  cases,  these  large  motions  are 

w 

mostly  waves  that  are  also  likely  to  be  controlled  by  the  local  value  of  N  if 

the  vertical  variation  of  N  is  small.   If  N(z)  varies  sharply,  then  other 

modes  can  be  generated  as  in  the  case  studied  by  Einaudi  and  Finnigan  (1981) . 

The  upper-air  soundings  in  Appendix  A  are  rather  different  to  their  situation. 

These  are  the  immediate  reasons  for  suspecting  that  0  ,  N,  and  d6/dz  are  the 

w 

main  parameters  determining  0„  and  the  heat  flux.   For  the  idealized  situation 
of  homogeneous  turbulence  in  a  strongly  stable  stratification,  the  Lagrangian 
model  of  Pearson  et  al.  (1982)  suggests  that 

°e  s  ce  (VN)(d^/dz)  »  (1'1) 

where  Co  is  a  temperature  fluctuation  parameter  of  the  order  of  unity.   In  fact 

u 

for  a  wide  class  of  turbulence,  £a  -  I/t/2.      The  model  also  predicts  that  the 
coefficient  of  thermal  diffusivity  is  given  by 

Kfl  --~L-«  F  <£/N   ,  (1.2) 

de/dz    w 

where  F  is  a  thermal  diffusivity  parameter  much  smaller  than  unity  and  roughly 

independent  of  G  / (L   N) .   The  relationship  in  (1.2)  can  also  be  written  to 
w   x 

show  that  the  effective  length  scale  of  turbulence  governing  heat  flux  (i.e., 
Kq/<j  )  is  equal  to  F  0  /N,   It  is  important  to  underline  that  this  theory 
makes  no  distinction  between  wavy  and  turbulent  motions  at  low  frequency,  and 
therefore  suggests  that  in  measuring  £fl  and  F  we  should  not  expect  any  great 
sensitivity  to  waviness. 

Consequently  in  Fig.  1.11  we  have  plotted  afl/(d0/dz)  as  a  function  of 

a  /N  for  all  three  of  the  days  and  conditions  considered  here  from  16  m  (where 

w  J 

d6/dz  is  evaluated)  to  275  m  (0„  is  interpolated).   We  have  only  included  points 
where  d6/dz  can  be  calculated.   In  some  cases,  d0/dz  is  indeterminate  because 
instruments  were  not  working  or  because  its  value  is  too  small  (0.0005°C/m) . 
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Figure  1.11. — Plot  of  the  relationship  in  eq.  (1.1)  showing  CQ 
respective  of  variations  in  z  or  L  or  the  existence  of  waves. 


-  0.8,  ir- 


In  the  latter  case,  temperature  fluctuations  persist,  but  there  can  no  longer 

be  a  local  connection  between  oQ  and  d0/dz.   (There  are  only  about  three  ex- 

w 

eluded  points,  and  all  are  at  225  m.) 

The  agreement  with  (1.1)  is  reasonably  good,  the  observed  constant  L,„ 
being  about  0.8  (see  Fig.  1.11),  compared  with  0.7  in  the  theoretical  model. 
Note  that  the  correlation  of  the  Kansas  surface  layer  data  for  strong  stabil- 
ity (z  <  25  m)  gives  £fi  between  0.5  and  1.0  (see  Appendix  B).   There  is  no 
apparent  trend  in  the  value  of  £q  with  z  or  with  the  existence  or  absence  of 
wave  motions. 


Figure  1.12a  shows  how  the  effective  length  scale  of  the  heat  flux  (Kq/°  ) 
varies  with  height  in  the  stable  boundary  layer.   Note  the  marked  variation  of 

the  reciprocal  heat  flux  a  /KQ  even  over  an  hour.   Figure  1.12b  shows  that  a  /K, 

w  6  &  w  t 

is  an  approximately  linear  function  of  N/a  ;  it  does  not  correlate  well  with 

the  velocity  gradient  length  scale  (dU/dz)/a  .   Finally  in  Fig.  1.13,  F  (i.e., 

w 

Kq/o  v  O  /N)  is  plotted  as  a  function  of  the  effectively  largest  scale  of 

motion  in  these  stratified  flows  (a  /N)  and  as  a  function  of  the  inverse 

W   -1/2 
square  root  of  the  Richardson  number  Ri     =  (dU/dz)/N.   There  is  no  obvious 

trend  in  this  data.   [According  to  surface  boundary  similarity  arguments  F 
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(a)  height  and  (b)  reciprocal  length  scales  imposed  by  buoyancy  N/a  and 
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w 


should  be  a  function  of  z/L  or  z/(G  /N) ,  when  z/L  >>  1/5.]   The  fact  that  the 

w 

average  value  of  F  is  about  0.17,  which  is  half  the  value  in  the  Kansas  ex- 
periment, might  be  explained  by  the  relatively  larger  values  of  a  associated 

w 

with  the  wave  motions  present  here.  However  there  is  no  direct  correlation 
between  wave  motion  and  any  of  the  results  in  Fig.  1.13. 

Recent  measurements  by  Nieuwstadt  (1982)  at  the  Cabauw  tower  in  the 
Netherlands  also  showed  that  in  stable  conditions  above  about  100  m  the  tem- 
perature fluctuations  and  heat  flux  could  be  scaled  on  the  local  turbulence 
and  temperature  gradient.   His  results  yield  £n  -  1.3  and  F  -  0.25.   Unlike  us, 
he  found  that  F  systematically  decreased  as  z  decreased,  reaching  a  value  of 
about  0.1  at  z/L  -  0.1. 
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1.5   TENTATIVE  CONCLUSIONS 

(1)   The  presence  of  coherent  wave  motions  in  stable  atmospheric  flows 
cannot  be  detected  from  the  mean  temperature  and  velocity  profiles,  nor  even 
on  some  occasions  from  spectra  of  the  vertical  velocity. 


(2)   Vertical  wave  motions  are  most  likely  to  be  detected  by  temperature 
spectra  and/or  temperature-vertical  velocity  cospectra.   Cospectra,  quad- 
spectra,  and  autocorrelations  help  identify  the  frequencies  where  wave  energy 
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is  particularly  concentrated.   Waves  of  higher  frequency  appear  to  be  cor- 
related over  smaller  vertical  distances.   Temperature  fluctuations  are  better 
correlated  over  larger  vertical  distances  than  are  the  velocity  fluctuations. 

(3)  It  may  be  useful  to  consider  turbulence  in  stably  stratified  condi- 
tions as  belonging  to  one  of  three  categories:  weak,  moderate,  or  strong  wave 
conditions,  depending  on  whether  significant  wave  effects  are  present  in 
neither  velocity  nor  temperature  spectra,  only  temperature  (or  w9  cospectra) , 
or  both  w  and  0  spectra,  respectively. 

(4)  Wave  motions  and  large-scale  turbulence  are  both  agents  for  the 
transport  of  heat  and  the  production  of  temperature  fluctuations.   In  the 
presence  of  a  stable  density  gradient,  where  the  buoyancy  frequency  N  is 
approximately  uniform  with  height,  it  appears  that  both  processes  are  asso- 
ciated with  vertical  motions  on  a  scale  approximately  equal  to  a  /N,  over  a 

-1  w 

time  scale  of  order  N   .   Consequently  dimensionless  parameters  can  be  defined 

to  estimate  (or  correlate)  temperature  fluctuations  G~  or  thermal  diffusivity 

KQ: 

tf0  ~  Ce  (Ow/N)(d6/dz) 

t)        WW 


where  typically 


Cn  -  0.8  +  0.25 
F  -  0.2  +  0.1  . 


[Note  that  these  general  ideas  are  not  appropriate  when  vertical  oscillations 
are  driven  by  elevated  inversions,  i.e.,  strong  variations  in  N.   Indeed,  as 
Finnigan  (CSIRO,  Canberra,  Australia;  personal  communication)  has  pointed  out, 
wave- induced  heat  flux  may  then  be  in  the  opposite  direction  to  the  high  fre- 
quency wave  and  turbulent  flux. ] 

(5)   The  integral  scale  for  the  vertical  turbulence  velocity  is  shown  to 
be  approximately  given  by 

rT(w),-l  M  dU/dz 

1  x  J   "a 
w 
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in  the  lower  part  of  the  boundary  layer,  at  z  -  25  m.   Below  this  height,  the 
direct  effect  of  the  ground  on  the  eddies  has  to  be  considered;  other  data  and 
theoretical  arguments  given  in  Appendix  D  suggest,  in  stable  neutral  and 
unstable  conditions, 

(w)  riB0-7duZd»  +  i 

x  a    z 

w 

In  the  upper  part  of  the  boundary  layer  and  above  it  (to  within  a  factor  of 
about  2) 

[L(w)]_1  -   N/0   (when  N  »  dU/dz)  . 
x  w 

(6)  The  integral  scales  of  temperature  L    and  w9  fluctuations  L     are 

(w)  X  X 

greater  than  those  of  L    .   (They  can  only  be  defined  in  weak  wave  conditions.) 

Typically  when  z  >  100  m  in  such  conditions  L(9)/L(w)  *  3;   l(w9)/L(w)  -  2. 

(w)  (6)   X    X         xx 

Near  the  surface,  L    decreases  where  L    does  not;  so  at  22  m  or  50  m, 

L<6>/L<">  .  7  +  1.  X 
xx       — 

(7)  whether  waves  are  present  or  not,  inertial  subrange  spectra  are 
observed  for  w  and  0.  Measurements  of  £   from  this  range  are  consistent  to 


within  a  factor  of  2  with  the  usual  estimate  based  on  the  turbulent  component 

2 
of  the  total  variance  (a  )  and  with  local  integral  scales  estimated  on  the 

basis  of  a  /N  and  the  wavenumber  at  the  peak  in  the  spectra. 

(8)  From  a  practical  point  of  view,  the  data  show  unmistakably  how  heat 
is  transferred  vertically  by  wave  or  wavelike  motions  as  well  as  by  low- 
frequency  turbulent  motions.   Consequently  contaminants  can  also  be  trans- 
ferred by  such  processes.   Therefore  it  is  important  that,  in  the  presentation 
and  recording  of  measurements  of  atmospheric  turbulence,  the  wavelike  motions 
should  be  included.   [By  excluding  spectral  contributions  from  frequencies 
below  0.001  Hz  (Caughey  et  al. ,  1979)  in  the  analysis  of  the  Minnesota  data, 
the  data  appeared  to  be  very  consistent,  presumably  because  it  was  just  a 
turbulence  record.   But  the  analysis  may  be  misleading  if  it  is  to  be  used  to 
estimate  heat  or  contaminant  fluxes.] 

(9)  The  second  practical  point  about  the  data  is  that  they  provide 
evidence  for  the  theoretical  suggestion  of  Pearson  et  al.  (1982)  that  in  a 
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stably  stratified  flow,  diffusion  heat  or  pollutant  transfer  takes  place  more 

by  slow  mixing  between  fluid  elements,  which  are  constrained  in  their  vertical 

motions  to  a  distance  of  about  O  /N,  rather  than  by  large-scale  advection  or 

w 

transport  by  fluid  elements.   If  F  -  0.2,  it  means  that  the  time  scale  for  a 
fluid  element  to  change  its  temperature  is  about  5  N  .   Therefore  vertical 
diffusion  of  any  pollution  released  from  an  elevated  source  into  the  atmo- 
sphere is  limited  to  a  /N  for  some  distance  downwind  of  a  source  until  mixing 
can  occur.   This  prediction  for  the  vertical  extent  of  a  plume  in  strongly 


stable  flows  can  also  be  estimated  from  the  variance  of  temperature  fluctua- 


2 
tions  afi  using  eq.  (1.1). 
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APPENDIX  A:  BACKGROUND  DATA  FOR  OBSERVATION  PERIODS 

Presented  here  are  background  data  on  synoptic  and  boundary  layer  scales 
for  the  three  nights  discussed  in  this  paper.   The  daily  weather  maps  for 
12  Z  (0500  MST),  published  by  NOAA's  Environmental  Data  and  Information  Ser- 
vice, are  reproduced  in  Figs.  A.1-A.3.   Two  maps  bracketing  the  selected 
nights  are  shown.   The  rawinsonde  plots  for  Denver  (30  km  from  the  BAO)  cover- 
ing the  same  observation  periods  are  presented  in  Figs.  A.4-A.6.   The  upper- 
air  soundings  at  1700  MST  the  evening  before  and  0500  MST  the  morning  after 
are  shown  to  highlight  changes  occurring  during  the  night.   On  a  smaller 
scale,  the  profiles  of  the  mean  and  turbulent  properties  of  the  stable  bound- 
ary layer  to  a  height  of  300  m  over  20-min  averaging  periods  can  be  found  In 
Tables  A.1-A.4. 
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Figure  A.l. — National  Weather  Service  maps  for  18  and  19  April  1978, 
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FRIDAY.  APRIL  21,  1978 


SATURDAY,  APRIL  22.  1978 


Figure  A. 2. — National  Weather  Service  maps  for  21  and  22  April  1978, 
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TUESDAY,  APRIL  15.  198(1 


WEDNESDAY.  APRIL  16.  1980 


Figure  A. 3. — National  Weather  Service  maps  for  15  and  16  April  1980, 
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APPENDIX  B:  VARIANCES  AND  SCALES  OF  TEMPERATURE  AND  VELOCITY 
FLUCTUATIONS  IN  THE  STABLE  BOUNDARY  LAYER* 

Monin-Obukhov  similarity  reasoning  and  measurements  over  flat  terrain 
show  that  the  mean  temperature  gradient  in  stable  conditions  (z/L  >  0)  is 
given  by 

T* 
d0/dz  =  —  (1  +  3  z/L),  (B.l) 

where  T^  =  -  w0/u^,  L  =  -u^  T  / (g  w0  k) ,  and  k  and  3  are  constants  (Monin  and 
Yaglom,  1971).   Since  fluxes  are  taken  as  constant  in  the  surface  layer, 
O   (-1.3  u^)  is  approximately  constant  with  height,  as  is  w0.   T  is  the  mean 
surface  temperature.   Consequently  for  similarity, 

re  =  ae/T*  =  constant,  when  z/L  >  0  .  (B.2) 

The  Kansas  data  give  a  range  for  Ffl  -  1.5  to  2.0  (Wyngaard  et  al.,  1971), 

whereas  Okamoto  and  Webb  (1970)  suggest  2.5  to  3.3,  whence  T  Q  =  -w0/GQa  -  0.35 

wU         o  w 

(Kansas  and  Minnesota  data,  Caughey  et  al.,  1979).   The  ranges  of  values  of 

TQ  and  T  .  measured  at  BAO  (z  <  22  m)  are  2.9  +0.6  and  0.23  +  0.1. 

0      wo  —  —  — 

Given  (B.l),  we  can  write  (B.2)  as  a  relation  between  the  rms  temperature 

—  (  Q  ^ 

fluctuation  G„,  the  mean  temperature  gradient  d0/dz,  and  a  length  scale  L 

defined  by 

00  =  L(9)  d0/dz  ,  (B.3a) 


where 


L<9>    -  i  +  1  z/l  re  ■  <B-3b> 


*  Summary  of  a  yet-to-be-published  note  by  J.C.R.  Hunt.  It  is  on  the  basis 
of  the  results  presented  here  that  some  of  the  measurements  are  discussed  in 
the  paper.   [Some  aspects  were  published  by  Hunt  (1982a).] 
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So  for  weak  stratification,  L    =  Tfi  kz,  and  for  strong  stratification,  where 
3  z/L  »  1, 


=   g   ~  (0.15  +  0.5)  L  , 


(B.3c) 


if  3  -  5.4,  as  in  the  Kansas  experiments  (Busch,  1973).   The  Monin-Obukhov 
length  scale  L  can  be  expressed  as 

r  3  T  _  -.1/2  r        J.    ,      ,„  x2  0  .,1/2 

L  = 


U*  L  To1 

_g(-w6)kj 


°l  (u*/°)2  £ 

w   "  w 


_(g/u^)(-w6/To)(3/kL)k  _ 


w 


(   31/2  ) 
\0.4  x  i.3y 


[1  +  l/(3z/L)] 


1/2 


(B.4a) 


Thus  for  strong  stratification, 


L  -  5  a  /N   for  the  Kansas  value  of  3. 
w 


(B.4b) 


At  BAO,  L/(a  /N)  varies  over  the  range  3  to  10  for  z  <  50  on  the  strongly 

w  1  - 

stratified  days  considered  here  when  z/L  >  — .   Thus  (B.3a)  becomes  for  the 
Kansas  data 


Oq    ~    V°w/N)  d6/dZ 


(B.5) 


where  Cfi  varies  from  0.5  to  1.0. 

2 
The  thermal  diffusivity  parameter  F  =  Kfl/(a  /N)  can  also  be  derived 

for  surface-layer  stable  turbulence.   From  (B.l),  the  coefficient  of  thermal 

diffusivity  is 


K, 


— Q  ku.   Z 

W0  * 


ku^  L 

■k 


d6/dz 


1  +  3  z/L   3[1  +  1/(3  z/L)]  * 


Then  from  (B.4a) 


K, 


F  = 


k(uA/0 ) 


w    3 


1/2 


a2/N   3[1  +  1/(3  z/L)]   N  k(a  /u*) 
w 


(1  +  1/3  z/L) 


1/2 


(1+3  z/L) 


1/2 
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g-1/2 
F  = 2   M        77Z  ,   /TN1l/2     *  (B'6a) 


w 

Thence 


(0„/u*)      [1  +  1/(3   (z/L)] 


F  -  — fTT  (B.6b) 

[1  +    (0/N)/«]1/Z 
w 


when  3(z/L)   »  1,   using  Kansas  data. 
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APPENDIX  C:  SCALES  OF  TEMPERATURE  AND  VELOCITY  FLUCTUATIONS 

If  in  stable  conditions  in  the  surface  layer  there  is  a  local  balance  of 

turbulent  kinetic  energy  (TKE)  between  production,  dissipation,  and  the  local 

flux  of  buoyancy,  then  the  steady-state  TKE  equation  becomes 

3  3 

u*      ?  u* 

£  =  i =  4  (dU/dz)  -  ££   ,  (C.l) 


£,U 


* 


which  provides  a  definition  for  the  energy  dissipation  Jcale  L    .   Since  it 

e'u* 
is  observed  that  in  stable  conditions 


dU/dz  -  T--  (1  +  az/L)  (C2) 


where  a  -  5  (Monin  and  Yaglom,  1971),  then  from  (CI)  and  (C.2) 


e  -  r-  [1  +  (a-Dz/L]  (C3a) 


(Busch,  1973),  and 


L     =  -  -    .   kj",  ._    .  (C.3b) 

e,uA   1  +  (a-l)z/L 

2 
The  implication  of  (C3a)  is  that  for  given  shear  stress  uA,  or  for  given 

2 
w  variance  a  ,  the  increase  of  dissipation  with  stable  stratification  is 
w  r 

(a  -  1),  or  about  4  times  greater  than  the  damping  of  turbulence  by  the 
buoyancy  flux.   Only  if  one  assumes  that  the  mean  horizontal  velocity  U  or 
mean  velocity  gradient  dU/dz  is  held  constant  as  z/L  increases  do  the  dissipa- 
tion and  turbulence  decrease.   Most  physical  discussions  of  the  damping 
action  of  stratification  are  rather  vague  as  to  what  aspect  of  the  flow  is 
kept  constant  or  determined  by  boundary  conditions  when  z/L  or  Ri  increases. 
An  explanation  for  this  result  is  given  in  Appendix  D.   [Usually  Prandtl's 
(1952,  p.  381)  discussion  of  "balls"  of  fluid  moving  vertically  influences 
most  qualitative  descriptions,  though  Turner  (1973,  p.  136)  emphasizes  that 
dissipative  viscous  losses  greatly  exceed  those  caused  by  buoyancy.   However, 
he  fails  to  consider  how  this  dissipation  changes  with  stratification  or  why 
it  does.]   Clearly  such  an  assumption  has  to  be  made  (implicitly  or  expli- 
citly) and  this  then  determines  which  is  the  dominant  physical  process. 
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Observations  and  physical  arguments  based  on  the  local  nature  of  small- 
scale  turbulent  motion  indicate  that 

e   =  C  a3/L(w) 
e  w  x 

even  in  stably  stratified  boundary  layers  where  C  is  a  coefficient  that 
varies  little  between  different  shear  flows  (see  Appendix  D) .   It  follows 
from  (C.l)  that 

l!w)  -  m°>*>3  \  „    .  <c-4a) 


x  e    w    *        e,u. 


whence   from   (C.3b) 


L(w)    «  C    (0   /uj\z/[l  +   (a-l)z/L]      .  (C.4b) 

x  £     w     * 


If   C^  =  0.4,    a,/u^  -  1.3,    and  k  =  0.4, 


w 


L(w)    =  0.35   z/(l  +  4z/L)    .  (C5) 

x 

The  results  in  (C.3a)  and  (C.5),  which  derive  from  the  observed  form  of 

<j>  -   kz(dU/dz)/u.  and  the  TKE  equation,  agree  well  over  the  range  of 
m  * 

0  <  z/L  <  1.0  with  Kaimal's  (1973)  Kansas  data.   (His  points  at  z/L  =  1.5  and 

(w) 
3.0  do  not  agree!)   This  equivalence  between  L     and  L   '  was  an  important 

£  >u&      x 
hypothesis  in  Townsend's  (1958)  argument  for  the  critical  Richardson  number; 

it  certainly  seems  to  be  substantiated  by  Kaimal's  data. 


The  relative  scale  of  temperature  fluctuations"  can  also  be  inferred  from 

vation; 
1958);  viz, 


2 
observations  and  the  equation  for  temperature  variance  0   (see  Townsend, 


-we  -^  +  f-  (we2/2)  -  efl    =  o  .  (c.6) 

dz   dz  t) 

2 
Then  let  £G  =  0Q  u./L   ,  where  L   is  the  temperature  dissipation  scale.   As 

9         ^  d   ~2 

in  the  TKE  equation  (C.l),  the  diffusion  term  ^—  (w0  /2)  is  negligible  in  the 

dz 

stable  surface  layer.   Since  d6/dz  is  observed  to  have  the  form  of  (B.l)  in 
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Appendix  B,  it  follows  that 

Qe  u*  °w   kZ 

\   =    (w^)2   (1  +  ^  Z/L)  '  (C'7) 

From  (C.3b)  and  (C.7)  the  ratio  of  temperature  to  vertical  velocity  dis- 
sipation scales  can  be  estimated  as 

i  /i    ..   l   /u*\l  +  (g-l)  z/l  ,  .. 

W0 

In  the  Kansas  data  T   Q  -  0.35,  and  a  /u.  -  1.3.   Assuming  that  the  tem- 
perature  integral  and  dissipation  scales  are  related  to  each  other  as  are  the 
vertical  scales,  then  when  z/L  =  0, 

L(6)  h 

x       * 
and  when  3  z/L  »  1,  since  3-6.4  (Busch,  1973), 


L(6)/L(w)  H  4>Q 
x    x 

In  the  Kansas  measurements  L   /L    -  8  when  z/L  =  0  and  3  when  z/L  -  1.0, 
whereas  in  the  BAO  measurements  (Kaimal  et  al.,  1982),  L   /L    -  5  when 
z/L  =  0.   These  results  give  some  support  to  the  arguments  leading  to  (C.9a) 

and  (C.9b).  We  conclude  that  there  may  be  at  most  a  reduction  by  a  factor  of  2 

(6)      (w) 
in  the  relative  scales  of  L    and  L    as  the  stability  increases  or  as  z 

xx 

increases  in  stable  conditions. 
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APPENDIX  D:  HYPOTHESES  CONCERNING  INTEGRAL  SCALES  AND 
DISSIPATION  SCALES  IN  SLOWLY  VARYING  TURBULENT  FLOWS 


The  ideas  in  this  Appendix  are  derived  from  studies  of  grid  turbulence 
(near  and  far  from  the  wall),  and  of  convective  and  neutral  atmospheric  bound- 
ary layers.   Some  of  the  new  results  are  applicable  in  these  other  flows  as 
well  as  in  stable  boundary  layers. 


Hypothesis  1 

In  many  turbulent  flows  it  is  observed  that 

£  «  c  a3/L(w)  .  (D.l) 

£  w  x 

In  turbulent  boundary  layers,  whether  stable,  neutral,  or  convective,  C 
varies  from  0.4  to  0.6,  whereas  in  grid  turbulence  C  -  1.5.   [For  some  justi- 
fication see  the  subsequent  analysis  as  well  as  Appendix  C,  and  Hunt  (1982b), 
Kaimal  (1973),  and  Townsend  (1958).] 


Hypothesis  2 


In  neutral  and  stable  shear  flows 


Ow  -  a_u^   ,   where   a  -  1.3  (D.2) 


(Townsend,  1976). 


(w) 
The  integral  scale  of  L  'is  defined  in  an  unconfined  stably  stratified 

shear  flow  by  dU/dz,  N,  and  a  .   A  linearized  analysis  similar  to  that  of 


w 


Townsend  (1976,  pp.  46-49)  for  the  rate  of  growth  T~  of  a  disturbance  with 
wavenumber  k  in  such  a  flow  shows  that,  if  |N/(dU/dz)|  «  1, 


+  VY"' 


T   =  X-tQOldU/dz  +  V(dU/dzr  -  X2(k)  N  |  ,  (D.3) 

where  the  constants  X-i >  X?  are  functions  of  k  and  hereafter  taken  to  be  equal 
to  1.0. 
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Weak  stable  stratification  affects  turbulent  energy  only  when  combined 
with  shear;  without  shear,  laboratory  experiments  of  decaying  turbulence  show 
(Britter  et  al. ,  1983)  that  there  is  little  effect  until  the  turbulence  is 
very  weak. 

If  the  shear  is  so  small  that  |dU/dz|  «  N,  and  the  turbulence  intensity 
is  weak  enough  that  NL   /a  >  3,  then  the  density  gradient  may  directly 

affect  the  straining  of  small  eddies  by  the  large  eddies,  and  thence  determine 

(w) 
£  and  L   .   Such  a  situation  is  likely  to  occur  only  above  the  boundary  layer 

X 

in  decaying  turbulence  or  in  turbulence  maintained  by  large-scale  wave  motion. 

For  these  reasons,  in  the  boundary  layer  where  dU/dz  -  N  and  NL   /a  <  1  there 

appears  to  be  no  physical  reason  why  the  turbulence  scale  should  be  controlled 

by  the  local  value  of  N/a  .   However,  when  the  flux  and  gradient  Richardson 

w 

numbers  are  constant,  as  assumed  by  Brost  and  Wyngaard  (1978)  in  their  model 

for  the  stable  boundary  layer,  dU/dz  -  N  and  then  the  local  value  of  a  /N  is 

(w) 
an  appropriate  scaling  for  L 

X 

2        2 
If  the  local  Richardson  number  Ri  =  N  /(dU/dz)   <  1/4,  then  it  follows 

from  (D.3)  that  the  natural  local  length  scale  L    in  an  unconfined  shear 

flow  is  given  by 

fwl  -1   t"1   ai  du/dz  /    X,  n2   \ 
w        w    \    4(dU/dz)  / 

where  A  is  an  unknown  constant. 

If  the  turbulence  is  close  to  a  rigid  boundary  at  z  =  0,  the  turbulence 
scale  is  reduced  in  proportion  to  z  (Hunt  and  Graham,  1978).   If  the  turbu- 
lence is  far  from  a  boundary  and  controlled  by  the  overall  scale  of  the  flow 
(e.g.,  a  grid  mesh  scale  L  or  a  wake  depth  L  or  a  convective  boundary  layer 

scale  L  or  the  scale  of  some  large-scale  wave  motion),  then  L  is  the  rele- 
o  o 

vant  scale. 

Thus  combining  these  suggestions  with  (D.4)  leads  to  the  next  hypothesis 
for  boundary  layer  flows  where  Ri  <  1/4. 
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Hypothesis  3 

i    ^      i        i    Ac  dU/dz/       ..2    \   A 

where  A  and  A  are  constants  for  the  shear  and  "blocking"  effects. 

From  the  kinematic  theory  of  shear-free  turbulence  near  a  rigid  interface 
(Hunt,  1982;  Hunt  and  Graham,  1978),  it  can  be  shown  that  close  to  the  sur- 
face when 

z/L  ->  0  or  [z/L(w)  (z  -»•  »)  ]  ■>  0,  then  L(w)  -  1.7  z  ,        (D.6) 

the  factor  1.7  being  determined  by  the  experimental  constant  in  Kolmogorov's 

inertial  subrange  law.   Since  in  the  limit  z/L  -*■  0  the  expression  in  (D.5) 

forTL   ]    equals  A/z,  it  follows  that 
x  B 

A,,  «  1/1.7  =  0.6  .  (D.7) 

a 

The  putative  constant  A  in  (D.5)  can  be  estimated  by  either  of  two 
alternative  observations. 


(1)   The  turbulent  kinetic  equation  in  a  boundary  layer  and  (D.l),  (D.2), 

and  (D.5)  show  that,  where  z/L  <<  1, 

o 

u^       0  /A„  dU/dz   A 


e  =  - —  =  C  a 
kz    e 


Since  it  is  observed  that 


Ik     dU/dz   A  V 


dU   u* 

dz^kz"  (1  +  az/L>  ■  (D.9) 

where  a  -  5  (Monin  and  Yaglom,  1971),  then  combining  (D.2),  (D.8),  and  (D.9), 
when  z/L  >>  1,  gives 

-3 


{h-A)~-°- 


As  =  kaJ  —^   _  aI  ^  0.7  if  C£  ^  0.6  (D.lOa) 


-   1.2  if  C  =  0.4  .  (D.lOb) 
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(2)  Alternatively,  the  value  of  A  can  be  fixed  by  comparing  the  expres- 

(w) 
sion  for  L    in  (D.5)  when  z/L  «  1  with  the  observation  that  in  a  neutral 
x     >  \  o 

boundary  layer  LW  »  (0.4  +  0.1) z  (Hunt  and  Weber,  1979).   In  that  case, 

Ag  a  1.3  x  0.4  x  [(2.5  +  0.5)  -  0.6]  (D.ll) 

=  1.0  +  0.25  . 

Thus  either  method  gives  roughly  equivalent  values  of  A  .   On  the  other  hand, 
if  C  -  0.6,  and  A  is  fixed  by  the  first  alternative,  then  (D.5)  provides  a 
prediction  that  in  a  neutral  SBL 

L(w)/z  *  0.5  . 
x 

(w) 
The  physical  implication  is  that  in  a  neutral  boundary  layer  L    is 

really  determined  by  two  comparable  effects — the  local  shear  and  the  wall 

proximity  effect. 

A  second  important  prediction  can  be  made  by  substituting  the  value  of 
the  "constant"  A  into  the  turbulent  kinetic  energy  equation  (C.l)  for  con- 
ditions near  the  ground  when  there  is  significant  stable  stratification,  i.e., 
z/L  >  0.   Then  from  (C.l),  (D.l),  (D.5),  and  (D.9) 


2  dU   u*   _ 
u*  d7  "  kL  "  Ce  ° 


so 


tf^W**  N2     V   VI 

WL   aw     \    4(dU/dz)2/   Z  J  ' 


dU   u*        (1  +  Sz/D 


dz   £ 

kZ  1  +  (A0  a^C  Ik 


(->  *4tis?) 


u*  -\ 

«  —  [1  +  (5  -  A.  a^C  /4)z/L]  ,  (D.12) 

Kz  S   3  £ 


where 


-   1  "  AS  a3  Ce         -      1 

k  »  1  ■         ■    and    a  -  ^ r   .  (D.13) 

C  aJ  B  A_  k  C  a, 

e  w  &        £  3 

Note  that  when  3z/L  »  1,  z/L  »  N2/(dU/dz)2,  so  dU/dz  -  —  (1  +   az/L).   Since 
A-  is  chosen  to  agree  with  the  measurements  of  the  neutral  case,  k  -  k  =  0.4, 
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but  3  is  not  fixed.   However  by  these  arguments  a   is  found  to  be  3.2.   Experi- 
ments over  flat  ground  give  a  value  of  a  equal  to  about  5 ,  which  is  certainly 
of  the  same  order  of  magnitude. 

We  see  in  Appendix  B  that  when  a  stable  stratification  is  applied  to  a 
turbulent  boundary  layer,  the  increase  in  viscous  dissipation  is  a  factor 
(a-1)  times  as  great  as  the  loss  of  energy  caused  by  the  buoyancy  flux.   We 
can  now  see  that  this  must  be  so.   For  given  uA  or  a  ,  the  generation  of  a 
small  buoyancy  flux  6B  must  require  an  increased  energy  input  by  an  increased 
shear  6(dU/dz).   If  this  were  the  only  balance  to  be  made,  then  6(dU/dz)  = 

6B/u..   But  an  increase  of  6(dU/dz)  reduces  the  integral  scale  L   '  and  in- 

2  X 

creases  the  dissipation  by  about  0.7  uA  6(dU/dz).   Thus  only  30%  (by  our 

calculation)  of  the  increased  shear  can  be  used  to  balance  the  buoyancy  flux. 

(w)  -1 

Only  by  realizing  the  relative  sensitivity  of  L   'to  both  dU/dz  and  z   can 

(    \ 

this  point  be  appreciated.   If  L    is  assumed  to  be  entirely  controlled  by 

the  distance  from  the  surface,  then  an  alternative  explanation  would  have  to  be 

sought  for  the  observed  sensitivity  of  £  to  z/L.   Since  the  turbulent  struc- 

(w) 
ture  at  scales  less  than  L    is  not  measurably  affected  by  stable  stratifi- 

x 

cation,  there  is  no  evidence  for  any  other  explanation.   [The  same  effects  are 
found  in  free  shear  layers  where  about  70%  to  80%  of  the  energy  of  Kelvin- 
Helmholtz  billows  is  dissipated,  the  rest  providing  a  buoyancy  flux  (Sherman 
et  al. ,  1978);  perhaps  the  same  explanation  is  appropriate.] 

As  a  third  prediction,  the  hypothesis  (D.5)  and  the  observation  (D.9) 
imply  that  when  az/L  >>  1,  dU/dz  =  auA/ (0. 4L) ,  and  that 

[L(w)fl  B  07   d£/dz  +  1  (D14) 

w 
or 

r  (w)  -1  a  6  x  Q.7   i 

L  x   J     0.4  x  1.3  l 

or 

L(w)  *  0.12  L  a  0.55  a  /N  ,  (D.15) 

x  w 

in  approximate  agreement  with  the  observations  (C.5)  and  (B.4b). 

Thus,  given  the  more  general  turbulence  hypotheses  (D.l),  (D.2),  and  (D.5), 
we  have  found  that  in  strongly  stable  conditions  L    must  be  proportional  to 
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L  and  O   /N  near  the  surface,  a  result  otherwise  only  deducible  from  the  kind 
w 

of  Monin-Obukhov  similarity  and  dimensional  arguments  used  in  Appendix  B. 


In  boundary  layers  with  strong  stable  stratification,  Brost  and  Wyngaard 

(1978)  argued  that  on  physical  grounds  a  /N  and  z  are  the  two  scales  that 

w 

determine  L    ,  so  by  interpolating  from  the  outer  edge  of  the  boundary  layer 

£>u* 
to  the  layer  nearest  the  surface,  they  postulated  that 

[L    J"1  =  N/a  +  z"1  .  (D.16) 

This  result  is  consistent  with  the  arguments  and  results  presented  here  when 
the  stratification  is  strong  and  at  the  outer  edge  of  the  boundary  layer  (even 

when  the  Richardson  number  is  not  constant).   When  the  stratification  is  weak, 

(w) 
(i.e. ,  L   '  N/a  <  1)  we  are  suggesting  on  the  basis  of  theoretical  arguments, 

X        w 

and  surface- layer  and  BAO  observations,  that  (D.16)  is  not  the  appropriate 
scaling,  even  in  boundary  layers.   The  more  general  scaling  suggested  here  is 
appropriate. 
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2.   WAVE  AND  TURBULENCE  STRUCTURE  IN  A  DISTURBED  NOCTURNAL  INVERSION 


Lu  Nai-ping 
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Bei j  ing ,  China 

W.D.  Neff  and  J.C.  Kaimal 
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Boulder,  Colorado  80303 


ABSTRACT.   Acoustic  sounder  and  tower  data  obtained  at  the  Boulder 
Atmospheric  Observatory  (BAO)  are  used  to  examine  several  features 
of  the  wave  and  turbulence  structure  associated  with  a  disturbed 
nocturnal  inversion.   General  features,  including  mean  fields  and 
Richardson  number,  for  the  case  selected  for  this  study  are  pre- 
sented.  Spectral  analysis  of  the  tower  data  reveals  a  separation 
of  energy  into  wavelike  and  turbulent  fluctuations.   Analysis  of 
the  heat  flux,  however,  shows  upward  counter-gradient  fluxes  in  the 
vicinity  of  a  low-level  jet  and  near  the  top  of  the  inversion.   Co- 
spectral  analysis  shows  that  the  major  contribution  to  the  upward 
heat  flux  occurs  at  frequencies  that  would  normally  be  considered 
characteristic  of  waves.   In  some  cases,  the  upward  flux  is  asso- 
ciated with  a  phase  shift  between  vertical  velocity  w  and  fluctuat- 
ing temperature  6  different  from  the  quadrature  relation  that  would 
be  expected  of  internal  waves.   Time  series  analysis  reveals  that 
these  unexpected  positive  fluxes  occur  in  relatively  short  bursts. 
Analysis  of  time  series  of  0  and  w  in  other  cases,  as  well  as  in- 
spection of  acoustic  sounder  records,  shows  that  sometimes  such 
upward  fluxes  can  result  from  a  combination  of  wave  motion  and 
horizontal  temperature  advection.   In  this  case  the  advection  is 
associated  with  a  shallow  cold  front. 


2.1   INTRODUCTION 

Remote-sensing  devices  using  microwaves  or  acoustic  waves  have  often  been 
used  to  provide  visual  documentation  for  special  case  studies  of  waves  and 
turbulence  within  the  nocturnal  inversion  layer  (e.g.,  Gossard  et  al. ,  1970; 
Emmanuel  et  al. ,  1972;  Merrill,  1977;  Einaudi  and  Finnigan,  1981).   whereas 
many  of  the  details  of  such  events  are  well  understood,  the  data  on  which  such 
analyses  are  often  based  arise  from  idealized  circumstances  (such  as  situations 
with  nearly  monochromatic  internal  waves).   More  commonly,  the  nocturnal 
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inversion  as  seen  by  remote-sensing  devices  reveals  a  complex  mixture  of  waves, 
instabilities,  and  turbulence.   Often,  the  complexity  of  the  inversion  seems  to 
be  related  to  particular  synoptic  situations  as  noted  by  Zhou  Ming-yu  et  al. 
(1980)  and  Neff  (1980).   The  acoustic  sounder  and  the  sensors  on  the  300-m 
tower  at  the  Boulder  Atmospheric  Observatory  (BAO),  operating  almost  con- 
tinuously over  the  last  4  years,  provide  the  data  base  needed  to  relate  typical 
patterns  in  nocturnal  inversion  development  to  synoptic  and  mesoscale  events. 

One  event  that  occurs  frequently  in  winter  at  BAO  is  the  passage  of  a 
shallow  cold  front.   Using  data  from  BAO,  we  present  a  case  study  of  one  such 
event.   We  will  first  describe  the  meteorological  setting,  the  general  features 
of  the  low-level  front,  and  the  complexity  of  the  waves,  turbulence,  and  ad- 
vective  processes  that  are  present.   We  will  then  outline  the  spectral  charac- 
teristics of  the  turbulence  and  waves  followed  by  detailed  interpretation  of 
the  features  contributing  to  an  anomalous  counter-gradient  heat  flux  that 
appears  in  a  conventional  20-min  covariance  of  w  (vertical  velocity)  and  6 
(fluctuating  temperature). 

2.2   SITE  CHARACTERISTICS  AND  INSTRUMENTATION 

BAO  is  located  about  30  km  north  of  Denver  and  30  km  east  of  the  foothills 
of  the  Rocky  Mountains.   This  site,  on  the  high  plains  of  Colorado,  is  charac- 
terized by  gentle  slopes  (Kaimal  et  al. ,  1982);  at  the  time  of  the  present 
study,  the  ground  cover  consisted  of  areas  of  wheat  stubble  mixed  with  bare 
areas  planted  with  winter  wheat. 

The  tower  instrumentation  at  BAO  has  been  described  in  some  detail  by 
Kaimal  and  Gaynor  (1983).   Briefly,  the  tower  is  instrumented  at  eight  levels 
(10,  22,  50,  100,  150,  200,  250,  and  300  m)  with  three-axis  sonic  anemometers, 
propeller-vane  anemometers,  fast-response  platinum  wire  and  slow-response 
quartz  thermometers,  and  a  cooled-mirror  dewpoint  hygrometer.   A  network  of 
five  sensitive  microbarographs,  measuring  pressure  fluctuations  p,  is  centered 
about  the  tower  to  detect  the  phase  speed  and  direction  of  internal  waves  pro- 
pagating across  the  site.   Two  acoustic  sounders  were  used:  one  close  to  the 
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tower  with  a  340-m  vertical  range,  and  a  second  about  600  m  away  with  an 
850-m  vertical  range. 

Data  from  the  tower  and  the  microbarographs  are  transmitted  via  phone 
line  to  the  central  data  archiving  computer  system  in  Boulder.   Standard 
system  programs  are  used  to  compute  fluxes,  variances,  and  20-min  spectra  in 
real  time.   Post-processing  programs  are  available  to  compute  spectra  as  well 
as  cospectra  and  quadrature  spectra  over  80-  or  160-min  periods.   A  special 
"beamsteering"  program  (Kaimal  and  Gaynor,  1983)  is  used  to  calculate  wave 
phase  speeds  and  directions  from  the  microbarograph  array.   This  program  has 
provisions  for  filtering  time  series  of  tower  data  in  any  specified  spectral 
interval  with  varying  passbands  and  cutoff  frequencies. 

2.3  GENERAL  DESCRIPTION  OF  CASE  STUDY 

2.3.1  Meteorological  Situation 

The  events  under  consideration  here  occurred  on  the  morning  of  27  March 
1981.   Earlier  in  the  evening,  the  development  of  the  nocturnal  inversion  had 
followed  an  uneventful  course  in  response  to  surface  cooling.   Winds  aloft 
were  from  south-southwest  reflecting  the  presence  of  a  cutoff  low  centered 
over  southwestern  Nevada.   The  combined  circulation  of  the  associated  surface 
low-  and  a  surface  high-pressure  system  located  over  the  Great  Lakes  main- 
tained southeasterly  surface  flow  over  the  eastern  plains  of  Colorado.   A 
surface  low  over  central  Wyoming  meanwhile  formed  and  moved  in  a  southerly 
direction  along  the  Continental  Divide  and  by  0500  MST  on  27  March  was  located 
south  of  Denver.   The  0500  MST  rawinsonde  ascent  made  by  the  National  Weather 
Service  at  Denver  showed  a  surface-based  inversion  capped  by  a  weak  isothermal 
layer  that  extended  from  800  to  700  mb.   The  winds  were  light  and  variable 
from  the  surface  to  700  mb,  above  which  the  winds  abruptly  increased  to  18  m/s 
from  the  southwest.   While  the  relatively  coarse  rawinsonde  network  revealed 
little  complexity  in  the  flow,  acoustic  sounder  and  tower  data  at  BA0  showed 
an  abrupt  change  in  the  nocturnal  inversion  shortly  after  midnight.   As  shown 
in  the  next  section,  winds  at  the  top  of  the  tower  became  light  from  the  west 
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while  in  the  lower  half  of  the  tower  a  10-m/s  jet  developed,  blowing  from 
north-northeast.   Associated  with  this  wind  shift  were  a  drop  of  about  3°C  in 
temperature  and  a  sudden  decrease  in  surface  pressure  of  about  0.5  mb.   Mea- 
surement of  the  wave  and  turbulence  properties  associated  with  this  unsteady 
jet  will  be  of  primary  interest  here. 

2.3.2  Acoustic  Sounder  and  Mean  Flow  Data 

Figure  2.1  shows  the  nocturnal  inversion,  as  seen  by  the  acoustic  sounder, 
during  the  period  from  midnight  to  0800  MST.   Three  distinct  periods  can  be 
noted.   The  first  is  the  initial  transition  period  between  midnight  and 
0100  MST  when  a  north-northeast  jet  first  appears.   The  second  is  between  0330 
and  0400  MST  when  the  inversion  is  strongly  disturbed.   The  third  is  after 
0600  MST,  just  before  the  onset  of  convection,  when  the  inversion  height  is 
relatively  stationary. 

The  intensity  of  the  echoes  depicted  in  Fig.  2.1  is  a  function  of  the 
strength  of  the  temperature  fluctuations.   Explicitly,  the  scattering  cross 
section  is  proportional  to  (Neff,  1980) 


ae 

8z 


5/3 

(2.1) 


The  first  term  in  brackets  measures  the  strength  of  the  turbulence,  where  K^ 
is  the  eddy  diffusion  coefficient  for  heat,  K^  is  the  coefficient  for  momen- 
tum, and  Ri  is  the  gradient  Richardson  number.   With  a  given  mean  potential 
temperature  (0)  profile,  the  distribution  of  turbulence  is  a  significant 
factor  in  the  echo  patterns  seen  in  the  acoustic  backscatter.   Results  from 
Neff  (1980)  also  suggest  that  the  distribution  of  wind  shear,  in  turn,  is  the 
most  significant  factor  determining  the  distribution  of  turbulence.   These 
conclusions  are  illustrated  fairly  clearly  in  Figs.  2.2  and  2.3.   Figure  2.2 
shows  the  wind  and  temperature  fields  as  well  as  isopleths  of  Ri  calculated 
using  data  from  fixed  levels  of  the  tower.   The  most  notable  feature  of  these 
figures  is  the  presence  of  two  scattering  layers:  one  near  the  ground,  the 
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0100 
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Figure  2.2. — Isopleths  of  wind  speed,  temperature,  and  Ri  developed  from 
tower  fixed-level  data  for  0000-0540  MST,  27  March  1981.   (Wind  speeds 
are  in  meters  per  second,  and  temperatures  are  in  degrees  Celsius.) 
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Figure  2.3. — Tower  wind  and  temperature  profiles  averaged  over  5  min  frc 
0105  to  0110  MST  compared  with  the  acoustic  sounder  record  from  0100  to 
0130  MST  obtained  at  a  site  just  160  m  north  of  the  tower.   The  dark, 
horizontal  bands  correspond  to  reflections  from  the  tower. 


second  elevated  above  a  clear  echo  region  on  the  record.   Comparison  of 
Figs.  2.1  and  2.2  indicates  that  the  clear  region  is  characterized  by  large 
values  of  Ri.   This  can  also  be  seen  in  a  detailed  comparison  between  the 
sounder  record  and  the  tower  wind  and  temperature  profiles  as  shown  in 
Fig.  2.3.   Within  the  resolution  permitted  by  the  50-m  spacing  of  fixed  tower 
levels,  the  non-echo  region  is  closely  associated  with  the  maximum  of  the  low- 
level  jet  where  the  shear  is  the  least.   The  elevated  scattering  layer  appears 
to  be  connected  with  the  region  of  increased  shear  and  stability  above  the 
jet.   In  accord  with  Eq.  (2.1),  the  strongest  scattering  occurs  in  layers  of 
largest  potential  temperature  gradient. 
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2.3.3   Internal  Wave  Behavior 

Of  interest  throughout  the  event  shown  in  Fig.  2.1  are  the  oscillations 
with  periods  of  several  minutes  at  the  height  of  the  inversion  base  capping 
the  low- level  jet.   Analysis  of  the  microbarograph  data  from  the  array  sur- 
rounding the  tower  shows  waves  of  several  frequencies  propagating  from  a 
generally  southwesterly  direction  at  speeds  from  6  to  22  m/s. 

To  determine  the  characteristics  of  these  waves,  we  used  the  BAO  beam- 
steering  program  to  calculate  wave  phase  speed,  direction,  and  coherence  as  a 
function  of  frequency.   This  program  requires  an  80-min  time  series.   Because 
waves  of  a  particular  frequency  may  occur  only  over  a  short  period  of  time,  we 
incremented  the  starting  time  of  the  calculation  by  20  min.   We  found  that 
prior  to  0300  MST,  only  one  wave  component  dominated  for  each  80-min  calcula- 
tion; that  at  the  time  of  frontal  passage,  the  wave  phase  speed  increased 
markedly;  and  that  after  0300  MST,  many  waves  of  different  speeds  and  direc- 
tions were  present.   These  results  are  shown  graphically  in  Fig.  2.4,  where 
data  are  plotted  at  the  center  time  of  each  period  and  the  vertical  bars 
represent  the  range  of  significant  wave  speeds  and  directions  found.   This 
broadening  of  the  wave  spectrum  coincides  with  the  large  disturbance  in  the 
inversion  shortly  after  0300  MST  as  seen  in  Fig.  2.1. 

Since  the  direction  of  propagation  of  these  waves  is  opposed  or  perpen- 
dicular to  the  low-level  flow,  it  may  be  that  the  waves  originated  from  a 
critical  level  aloft.   The  most  likely  region,  determined  from  the  0500  MST 
rawinsonde  launched  from  Denver,  is  at  2  km  where,  in  a  strongly  sheared  re- 
gion, the  winds  increase  from  less  than  1  m/s  to  nearly  18  m/s.   This  seems  a 
reasonable  conclusion  based  on  the  work  of  Keliher  (1975)  who  studied  the 
correlation  of  surface-measured  waves  with  jet-stream  velocities  and  direc- 
tions.  Atkinson  (1981)  suggests,  though,  that  although  a  low-level  inversion 
may  amplify  surface  pressure  fluctuations,  wavelengths  less  than  5  km  at 
tropopause  level  are  unlikely  to  be  associated  with  surface  pressure  fluctua- 
tions.  Prior  to  0300  MST,  the  dominant  wave  period  was  about  15  min;  after- 
wards wave  periods  ranged  from  4  to  15  min.   On  the  basis  of  the  wave  phase 
speeds  in  Fig.  2.4  after  0300  MST,  only  the  long-period  oscillations  in  the 
inversion,  from  10  to  20  min  in  length,  are  likely  to  be  associated  with  an 
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Figure  2.4. — Wave  phase  speed  and 
direction  calculated  at  20-min  inter- 
vals from  1900  MST  on  26  March  to 
0800  MST  on  27  March  1981.   Calcula- 
tion is  in  overlapping  80-min  seg- 
ments.  Vertical  bars  represent  the 
range  of  significant  wave  speeds  and 
directions. 
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upper-level  source.   The  origin  of  the  shorter  period  waves  is  not  determined 
although  one  might  speculate  that  they  could  be  an  artifact  of  the  frontal 
disturbance  and/or  its  interaction  with  the  foothills  of  the  Rocky  Mountains 
some  30  km  to  the  west. 


2.4   SPECTRAL  ANALYSIS 

2.4.1  Spectra  Under  Steady  Conditions 

The  evidence  of  waves  in  both  the  microbarograph  and  acoustic  sounder 
data  prompted  a  closer  examination  of  the  tower  wind  and  temperature  data  by 
spectral  techniques.   The  first  period  chosen  for  this  analysis  was  one  with  a 
fairly  constant  inversion  height,  just  before  the  onset  of  convection. 
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A  number  of  authors  have  discussed  spectra  under  stable  conditions. 
Steward  (1969)  pointed  out  some  of  the  difficulties  of  distinguishing  rigor- 
ously between  waves  and  turbulence.   However,  a  number  of  subsequent  papers 
(Axford,  1971;  Kaimal,  1973;  Caughey  and  Reading,  1975;  Caughey,  1977;  Caughey 
et  al. ,  1979)  have  shown  that  in  the  analysis  of  power  spectra  there  is  often 
a  clear  separation  into  wavelike  and  turbulent  regions.   Figure  2.5  presents 

such  results  from  data  taken  between  0620  and  0720  MST,  27  March.   These  data 

_2 
show  a  clear  spectral  gap  at  about  10  "Hz.   In  the  w  spectra,  the  gap  is  much 

less  noticeable  near  the  ground,  primarily  because  w  must  approach  zero  at  the 

surface.   The  high-frequency  portion  of  the  spectra  follow  the  characteristic 

-2/3  power  law.   Cross-spectra  at  250  m  for  the  same  time  period  show  an 

-2 
abrupt  change  in  spectral  slope  at  10  'Hz.   The  slope  at  lower  frequencies  is 

almost  -5;  the  slope  at  higher  frequencies  tends  to  be  a  bit  less  than  -A/3 

predicted  from  surface  layer  similarity  theory  (Wyngaard  and  Cote,  1972). 

Integrating  the  wavelike  and  turbulent  portions  of  the  spectra  sepa- 
rately, we  get  the  results  in  Table  2.1.   Here  it  can  be  seen  that  the  low- 
frequency  portion  of  the  spectra  and  cospectra  contributes  more  strongly  to 
the  variances  and  covariances  than  the  high  frequency  portion  does. 

Table  2.1 — Ratios  of  low-frequency  (L)  to  high-frequency  turbulent  (H) 
contributions  to  standard  deviations  and  covariances  during  the 
period  (0620-0720  MST,  27  March  1981) 
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2.4.2  Disturbed  Cases 

Figure  2.6  shows  time  series  of  temperature  from  2100  to  0540  MST,  and 
Fig.  2.7  shows  time  series  of  temperature  flux  w0  from  0000  to  0540  MST. 
(Prior  to  0000  MST,  the  wind  was  through  the  tower  rendering  the  sonic  mea- 
surements useless.)   Data  for  these  time  series  were  obtained  as  20-min  means, 
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Figure  2.6. — Time  series  of  temperature  averaged  over  20  min  at  four 
tower  levels  from  2100  MST  on  26  March  through  0540  MST  on  27  March. 


The  covariance  of  temperature  and  vertical  velocity  was  computed  in  the  usual 
manner  after  removal  of  the  means.   The  temperature  and  temperature  flux  data 
used  are  the  standard  archived  values  in  the  BAO  files.   As  we  will  show 
below,  use  of  these  standard  flux  values  requires  some  caution.   In  principle, 
simple  removal  of  mean  values  should  be  satisfactory  if  no  more  than  one  time 
series  has  a  trend.   In  general  it  is  assumed  that  the  time  series  of  w  is 
free  of  trends.   An  example  will  be  given  later  where  this  is  not  the  case. 

Notable  in  the  temperature  traces  is  the  interruption  of  the  general 
cooling  trend  by  a  warming  at  midlevels  after  0300  MST.   Unusual  in  the  flux 
time  series  are  the  positive  fluxes  observed  several  times  during  the  night. 
Since  these  are  counter-gradient,  they  are  most  puzzling.   Applying  the  tech- 
nique used  in  the  previous  section,  of  integrating  the  cospectrum  of  w  and  9 
separately  over  wave  and  turbulence  intervals,  we  find  that  the  positive  con- 
tributions to  the  temperature  flux  come  from  low-frequency  motions.   Inspec- 
tion of  Fig.  2.7  shows  positive  fluxes  mostly  in  the  midlevels  of  the  tower 
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Figure  2.7. — Time  series  of  temperature  flux  w0  calculated  over  20  min 
from  0000  to  0540  MST,  27  March  1981. 


(100-200  m)  with  negative  fluxes  at  250  and  300  m.   Comparison  with  the 
temperature  traces  in  Fig.  2.6  shows  that  this  midlevel  flux  convergence 
corresponds  to  a  period  of  net  warming  at  200  m.   From  Fig.  2.7,  taking  an 
average  flux  difference  of  0.2°C  m/s  between  150  m  and  250  m  respectively, 
over  the  period  from  0320  to  0440  MST,  produces  a  net  heating  of  9.6°C  in  the 
absence  of  advective  or  radiative  temperature  changes.   This  temperature 
change  compares  with  the  observed  heating  of  about  4°C  at  200  m  over  the  same 
period.   In  the  next  section  we  examine  these  unusual  counter-gradient  fluxes 
in  more  detail. 
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2.5   INTERPRETATION 

Our  analysis  of  the  temperature  flux  measured  at  each  level  showed 
periods  during  which  there  was  a  net  upward  flux  that  was  counter  to  the  mean 
gradient.   The  first  of  these  occurred  shortly  after  0100  MST,  the  second 
after  0500,  and  the  third  after  0330  MST. 

In  the  analysis  of  the  w0  cospectrum,  upward  temperature  flux  was  as- 
sociated with  frequencies  characteristic  of  internal  waves.   The  spectral 
technique,  however,  has  the  limitation  that  it  does  not  allow  an  identifica- 
tion of  contributions  to  the  positive  flux  from  any  isolated  events.   In 
addition,  it  is  difficult  to  separate  idiosyncracies  in  the  time  series  (such 
as  trends  or  steps)  from  waves  after  they  have  been  Fourier- transformed. 

Inspection  of  the  time  series  of  the  product  w0  often  showed  singular 
flux  events  where  an  excursion  in  one  direction  was  followed  by  one  of  the 
opposite  sign.   Since  the  net  contribution  appeared  to  come  from  the  imbalance 
between  the  two,  which  was  not  always  evident  from  the  actual  time  series  of 
w6,  a  time  series  of  the  cumulative  sum  of  wO  is  more  useful  for  identifying 
the  contribution  of  such  individual  events.   Such  time  series,  plus  those  of  w 
and  6,  are  shown  for  each  of  three  cases  in  Fig.  2.8.   Case  1  shows  the  effect 
of  a  long-period  wave  combined  with  a  step  decrease  in  temperature  associated 
with  the  passage  of  the  frontal  interface.   Case  2  shows  a  similar  case  with 
the  added  complication  of  a  slow  drift  in  the  mean  vertical  velocity.   Case  3 
represents  a  strong  disturbance  in  the  inversion. 


2.5.1  Effects  at  the  Frontal  Interface,  0100-0140  MST 

A  key  feature  in  Case  1  is  the  sharp  step  in  the  contributions  to  the 
flux  associated  with  the  time  when  the  inversion  interface  moves  through  the 
150-m  tower  level.   As  indicated  by  the  vertical  velocity  trace,  the  passage 
occurs  at  the  particular  phase  of  the  wave  so  that  a  negative  velocity  is 
associated  with  a  decrease  in  temperature.   Were  the  interface  to  remain 
stationary  and  only  perturbed  by  a  wave,  there  would  be  no  net  flux.   However, 
with  advection  of  the  tilted  frontal  interface  and  its  associated  waves,  this 
becomes  a  singular  event. 
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2.5.2   Trends/Steps  in  the  Data,  0420-0540  MST 

Case  2  shows  a  similar  set  of  traces  for  the  period  after  0500  MST.  In 
this  case  there  is  a  gradual  increase  in  the  net  positive  flux  after  an  ini- 
tial step.  An  inspection  of  the  time  series  of  w,  however,  shows  a  downward 
trend  in  w,  of  about  0.25  m/s  over  40  min,  after  the  abrupt  drop  in  tempera- 
ture at  0450  MST.  This  trend,  together  with  an  offset  in  temperature,  con- 
tributes the  net  upward  flux.  Separate  calculations,  based  on  means  taken 
before  and  after  the  step,  show  no  positive  flux. 


2.5.3   Strong  Disturbances  in  the  Inversion,  0320-0440  MST 

Another  anomalous  case  (Case  3)  that  we  noted  earlier  occurred  between 
0320  and  0400  MST.   Again,  calculating  the  cumulative  contributions  to  the 
temperature  flux,  we  see  that  the  contributions  occur  in  discrete  steps 
(Fig.  2.8).   However,  inspection  of  the  time  series  shows  no  steps  or  trends 
in  either  the  temperature  or  vertical  velocity  data  at,  for  example,  150  m  as 
shown  in  Fig.  2.8.   Returning  to  Fig.  2.7,  it  should  also  be  noted  that  the 
fluxes  at  250  and  300  m  due  to  the  same  disturbance  are  downward. 

Spectral  analysis  of  w  and  0  from  150  m  reveals  several  peaks  in  the 
range  from  0.001  to  0.004  Hz  (Fig.  2.9).   Use  of  a  linear  frequency  scale  here 
allows  identification  of  three  distinct  frequency  components  common  to  both  w 
and  0.   Of  note  is  that  the  central  peaks  match  well  while  the  adjacent  peaks 
are  shifted  in  frequency.   The  cospectrum  of  w  and  0,  using  standard  BAO 
programs,  shows  a  broad  peak  in  the  same  region.   However,  logarithmic  smooth- 
ing of  the  spectra  by  these  programs  does  not  allow  detailed  examination  of 
the  cospectrum.   Averaging  of  the  cospectra  and  quadrature  spectra  over  the 
entire  band  gives  an  average  phase  shift  of  -62°. 

Not  shown  in  Fig.  2.9  is  the  spectrum  of  pressure  fluctuations  at  the 
surface.   In  this  case  the  pressure  spectrum  closely  followed  the  w  spectrum. 
From  our  earlier  analysis  of  microbarograph  data,  such  as  that  presented  in 
Fig.  2.4,  we  also  found  a  spectrum  of  waves  propagating  opposite  in  direction 
to  the  inversion  jet.   To  correlate  surface  pressure  fluctuations  with  events 
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Figure  2.9. — Spectra  of  w  and  0  from 
150  m  in  the  frequency  range  1-4  x 

-3 
10   Hz  on  linear  scales,  for  0320- 

0440  MST,  27  March  1981. 


Figure  2.10. — Cospectrum  of  w  and  p 
for  the  time  period  in  Fig.  2.9. 


within  the  inversion,  we  computed  the  cospectrum  of  w  at  150  m  with  surface 
pressure  at  the  base  of  the  tower.   This  is  shown  in  Fig.  2.10  where  smoothing 
of  the  spectrum  has  again  produced  one  central  peak  in  the  cospectrum.   The 
quadrature  spectrum  in  this  case  is  small  compared  with  the  cospectrum.   This 
in-phase  relationship  may  reflect,  in  part,  a  phase  shift  introduced  by  the 
slower  time  constant  of  the  pressure  sensor  compared  with  that  of  the  sonic- 
derived  w  measurement.   The  finding  that  variations  in  w  within  the  inversion 
are  well  correlated  with  surface  pressure  fluctuations  remains  valid.   Similar 
comparisons  with  temperature  are  not  as  clear  except  in  those  cases  where  the 
temperature  sensor  happened  to  be  near  the  top  of  the  elevated  inversion  (see 
Fig.  2.3)  where  temperature  fluctuations  are  more  pronounced. 

Given  the  presence  of  these  cospectral  peaks,  we  used  our  filtering 
program  to  produce  the  time  series  of  w  and  0  in  the  frequency  range  shown  in 
Fig.  2.9.   These  are  shown  in  Fig.  2.11  together  with  identifiable  phase 
shifts  at  each  cycle  of  the  wave  at  heights  of  300  and  200  m.   The  fact  that 
the  wave  packet  centered  on  the  time  of  the  primary  disturbance  is  fairly 
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Figure  2.11. — Narrow  bandpasses  time  series  of  w  and  G  at  two  fixed 
levels  for  0320-0440  MST  on  27  March  1981  for  the  frequency  range 
shown  in  Fig.  2.9.   <f>  indicates  individual  phase  shifts. 


broad  is  a  consequence  of  the  narrow  bandwidth  used  to  filter  the  original 
time  series.   The  breadth  of  the  packet,  At  in  time,  obeys  the  relation 
Af'At  =  1,  where  Af  is  the  bandwidth.   For  an  ordinary  internal  wave,  density 
variations  are  90°  out  of  phase  with  vertical  velocity  (e.g.,  Turner,  1973). 
As  can  be  seen  in  comparison  with  Fig.  2.7,  a  shift  of  less  than  90°  leads  to 
a  positive  flux,  whereas  a  shift  greater  than  90°  produces  a  negative  flux. 

Although  these  results  suggest  vertical  transfer  of  heat  at  wavelike 
time  scales,  the  actual  transfer  of  heat  by  a  wave  seems  physically  implausi- 
ble.  However,  one  can  imagine  a  combination  of  circumstances  involving  large- 
amplitude  solitary  waves  acting  in  conjunction  with  horizontal  temperature 
advection  to  produce  such  isolated  fluxes.   Profiles  at  a  single  location  are 
not,  however,  sufficient  to  resolve  such  questions. 
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2.6   CONCLUSIONS 

The  analysis  of  a  nocturnal  Inversion  disturbed  by  the  passage  of  a  shal- 
low cold  front  shows  the  dominance  of  low-frequency  motions  over  turbulence  in 
the  calculation  of  variances  and  covariances.   In  addition,  the  unsteadiness 
of  the  flow,  with  both  wave  and  advective  processes  present,  was  shown  to 
create  difficulties  in  standard  flux  calculations.   Such  effects  are  critical 
in  the  interpretation  of  long-term  flux  data  archived  from  towers  such  as  BAO. 

Three  cases  were  presented.   The  first  one  showed  that  anomalous  fluxes 
could  occur  with  the  fortuitous  superposition  of  wave  motion  and  the  advection 
of  the  frontal  interface  through  a  fixed  tower  level.   The  second  revealed  a 
positive  flux  arising  from  a  combination  of  a  trend  in  w  together  with  a  sharp 
decrease  in  temperature.   The  third,  involving  a  large-amplitude  disturbance 
in  the  inversion,  proved  more  intriguing.   Analysis  of  this  third  case  moti- 
vated the  introduction  of  two  additional  analysis  techniques.   The  first 
involved  forming  a  cumulative  sum  of  the  product  of  w  and  0.   This  showed  that 
the  apparent  wavelike  contributions  to  the  upward  vertical  temperature  flux 
actually  occurred  in  short  bursts.   Cospectral  analysis  of  w  and  6  in  this 
case  showed  equal  and  opposite  phase  shifts,  away  from  the  90°  expected  of 
internal  wave  motion,  above  and  below  the  elevated  inversion  base.   A  second 
technique  was  introduced  to  relate  events  within  the  inversion  to  wave  motion 
sensed  by  the  surface  microbarograph  array.   To  accomplish  this,  we  formed 
cospectra  of  tower  variables  such  as  w  with  the  pressure  time  series  at  the 
base  of  the  tower.   The  cospectrum  of  w  and  p  showed  strong  peaks  at  the  same 
frequencies  as  those  in  the  spectra  of  w  and  0,  suggesting  that  the  large- 
amplitude  disturbances  in  the  inversion  layer  flow  were  wave  related.   Al- 
though some  of  the  tower  data  show  that  advection  may  be  an  important  element 
in  the  interpretation  of  such  low-frequency  effects,  sufficient  data  to  re- 
solve the  problem  do  not  exist. 
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ABSTRACT.   The  structure  of  nocturnal  inversions  in  the  first  300  m 
of  the  atmosphere  is  analyzed  using  the  observational  data  from  the 
Boulder  Atmospheric  Observatory  from  March  through  June  1981.   The 
temperature  profiles  show  more  than  one  inversion  layer  41%  of  the 
time  during  the  observational  period.   The  vertical  distributions  of 
wind  speed  and  moisture  also  show  evidence  of  stratification  during 
these  multiple-layer  events.   The  relation  between  the  radiative 
cooling  rate  in  time  and  height,  including  moisture,  and  the  vertical 
structure  of  the  multiple  layers  is  calculated.   The  vertical  dis- 
tribution of  eddy  kinetic  energy  and  the  turbulent  vertical  fluxes 
of  heat  and  momentum  are  also  calculated.   Turbulent  structure  in 
the  elevated  inversion  layers  is  more  complicated  than  that  in  the 
single-layer,  stable  nocturnal  boundary  layer.   The  total  heat  bud- 
get for  a  multiple-layer  case  is  calculated,  and  turbulent  cooling 
is  found  to  be  negligible  relative  to  radiative  cooling  and  to  hori- 
zontal advection  and/or  horizontal  divergence  of  heat  flux. 


3.1   INTRODUCTION 

In  recent  years,  the  study  of  the  nocturnal  stable  boundary  layer  has 
received  much  attention.   Many  nocturnal  boundary  layer  models  have  been 
developed  (Delage,  1974;  Businger  and  Arya,  1974;  Wyngaard,  1975;  Blackadar, 
1976;  Brost  and  Wyngaard,  1978;  Andre  et  al.,  1978;  Zeman,  1979;  Li  et  al., 
1982a),  and  many  studies  have  been  made  of  the  parameterization  of  the  noc- 
turnal boundary  layer  as  well  as  of  the  basic  vertical  structure  of  the 
stably  stratified  atmospheric  boundary  layer  (Rao  and  Snodgrass,  1979; 
Sundararajan,  1979;  Mahrt  et  al.,  1979).   However,  observational  studies  of 
the  nocturnal  lower  atmosphere,  of  which  the  boundary  layer  is  only  a  part, 
are  rather  incomplete.   Its  description  has  been  restricted  because  of  the 
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lack  of  detailed  turbulence  and  mean  field  measurements  with  height.   The 
Boulder  Atmospheric  Observatory  (BAO)  tower  provides  such  information. 

An  experiment  directed  toward  the  study  of  the  nocturnal  lower  atmo- 
sphere was  conducted  at  BAO  between  March  and  June  1981.   Measurements  taken 
during  this  intense  observational  period  included,  in  addition  to  standard 
data  from  the  300-m  BAO  tower,  profiles  of  wind,  temperature,  and  humidity 
using  instruments  attached  to  a  movable  carriage  on  the  tower,  and  facsimile 
records  of  echo  intensities  from  an  acoustic  sounder  located  near  the  tower. 

The  sensors  on  the  tower  and  data  acquisition  procedures  have  been  dis- 
cussed by  Kaimal  and  Gaynor  (1983).   The  carriage  sensors  were  mounted  on  two 
horizontal  booms  separated  10.3  m  vertically.   The  lower  boom  supported  a 
sonic  anemometer  measuring  the  mean  and  fluctuating  components  of  the  wind 
along  three  orthogonal  axes,  a  quartz  thermometer  measuring  mean  temperature, 
a  platinum-wire  thermometer  measuring  temperature  fluctuations,  a  Lyman-alpha 
hygrometer  measuring  fluctuations  in  specific  humidity,  and  a  propeller  vane 
anemometer  measuring  mean  wind  speed  and  direction.   For  gradient  measure- 
ments across  the  10.3-m  vertical  spacing,  the  upper  boom  supported  another 
quartz  thermometer,  Lyman-alpha  hygrometer,  and  propeller-vane  anemometer. 
The  movable  carriage  ascends  the  full  tower  length  in  about  8.8  min  and 
descends  in  about  8.5  min. 


3.2   FORMATION  OF  MULTIPLE  LAYERS 

Measurements  made  with  the  acoustic  sounder  and  "with  sensors  mounted  on 
the  tower  indicate  that  the  structure  of  the  nocturnal  lower  atmosphere  is 
very  complex.   The  inversion  structure  is  often  correlated  with  isolated  echo 
layers  in  the  facsimile  records.   The  statistical  results  in  Fig.  3.1  show 
that  of  the  87  nights  examined  (25  March-22  June  1981),  35  had  at  least  one 
occurrence  of  multiple  layers.   The  probability  of  occurrence  of  multiple 
layers  on  any  one  night,  therefore  seems  reasonably  high.   The  starting  time 
of  such  events  is  variable,  but  Fig.  3.1  shows  a  high  concentration  (86%  of 
the  cases)  between  2000  and  0200  mountain  standard  time  (MST).   The  duration 
of  a  multiple-layer  event  varies  from  3  to  12  h. 
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Figure  3.1. — Frequency  distribution 
of  times  when  multiple  inversion 
layers  initially  occurred  between 
25  March  and  22  June  1981. 
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The  formation  of  multiple  layers  is  almost  always  associated  with  cer- 
tain atmospheric  conditions  and/or  local  topographical  features.  Following 
are  some  of  the  conditions  that  favor  such  formation: 

(1)  Subsidence  and  advection  associated  with  the  passage  of  a  cold 
front  in  uniform  and  complex  terrain  (Neff,  1980). 

(2)  Cold-air  outflows  from  local  thunderstorms  (Goff,  1976). 

(3)  Drainage  winds  in  the  vicinity  of  mountain  ranges,  as  observed 
often  at  BAO  (Hootman  and  Blumen,  1981). 

Li  et  al.  (1982b)  have  discussed  a  multiple-layer  event  caused  by  advec- 
tion.  This  study  will  focus  on  a  more  typical  event  at  BAO,  one  caused  by 
drainage  flow  from  the  Rocky  Mountains  60  km  to  the  west. 


3.3  A  CASE  STUDY 

The  event  under  consideration  occurred  between  0100  and  0530  MST, 
24  April  1981.   One  carriage  profile  (either  ascending  or  descending)  was 
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Figure  3.3. — The  time-height  cross  section  of  the  wind  direction  on 
24  April  1981.   Isolines  are  at  30°  intervals. 


recorded  every  10  min  during  that  period.   The  observed  time  series  on  all  the 
data  channels  were  block-averaged  over  consecutive  10-s  blocks.   This  averag- 
ing corresponds  to  a  6-m  spatial  averaging  in  the  vertical  and  is  clearly  more 
useful  for  detection  of  elevated  layers  than  the  50-m  spacing  between  the 
fixed  levels  on  the  tower.   These  smoothed  vertical  profiles  from  the  carriage 
sensors  were  used  to  construct  time-height  cross  sections  of  the  measured 
variables. 


The  acoustic  record  (Fig.  3.2)  indicates  that  multiple-layer  structures 
started  to  form  between  0100  and  0200  MST  and  between  0230  and  0330  MST.   The 
time-height  cross  section  of  wind  direction  between  0120  and  0530  MST  on 
24  April  1981  is  shown  in  Fig.  3.3.   Before  0230  MST,  the  wind  direction  was 
essentially  southerly  above  150  m,  northwesterly  between  150  and  50  m,  and 
quite  variable  below  50  m.   After  0230  MST,  a  significant  change  in  wind 
direction  occurred  within  the  tower  height.   The  wind  direction  became  gen- 
erally northwesterly  above  50  m  and  southeasterly  below  50  m. 
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0530 


Figure  3.4. — The  time-height  cross  section  of  the  temperature  on  24  April 

1981.   Solid  lines  are  isotherms  at  0.2°C  intervals;  thin  dashed  lines 

are  at  0.5°C  intervals;  thick  dashed  lines  indicate  tops  of  the  inversions, 


The  time-height  cross  section  of  temperature  for  the  same  case  is  shown 

3T 
in  Fig.  3.4.   The  thick  dashed  lines  indicate  tops  of  the  inversions  (tt—  >  0 

dz 

below  the  thick  dashed  lines).   The  temperature  field  shows  multiple  inver- 
sion layers  between  0120  and  0200  MST,  in  agreement  with  the  multiple  echo 
layers  on  the  acoustic  record  (Fig.  3.2).   After  0230  MST,  there  was  an 
inflow  of  cold  air  within  the  middle  and  lower  tower  levels.   The  inversion 
then  divided  into  two  parts,  the  lower  part  continuing  to  exist  as  a  radia- 
tion inversion  at  50  m,  the  upper  part  lifting  to  250  m  by  0300  MST.   The 
temperature  within  the  upper  inversion  increased  about  1°C  with  height  in  a 
20-  to  30-m  thickness.   These  layers  correspond  to  strong  echo  layers  in  the 
acoustic  record  of  Fig.  3.2.   The  thermal  stratification  between  the  two 
inversions  was  nearly  isothermal.   After  0300  MST  an  elevated  inversion  seems 
to  hover  around  150  and  200  m  as  seen  in  Fig.  3.4.   Although  Fig.  3.4  shows 
only  two  inversion  layers,  the  acoustic  sounder  record  in  Fig.  3.2  displays 
at  least  four  echo  layers  above  50  m.   The  multiple  echo  layers  may  represent 
the  top  and  bottom  of  each  inversion  layer,  or  they  may  represent  a  smaller 
scale  structure  in  the  temperature  profile  not  easily  resolvable  in  Fig.  3.4. 
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Figure  3.5. — The  time-height  cross  section  of  the  wind  speed  on  24  April 
1981.   Solid  lines  are  isotacks  at  0.5-m/s  intervals;  thick  dashed  lines 
indicate  maximum  values  of  wind  speed,  i.e.,  the  axes  of  the  low-level 
jets. 


The  time-height  cross  section  of  wind  speed  for  this  case  is  shown  in 
Fig.  3.5.   The  wind  field  shows  multiple  wind  speed  maxima  from  0120  to  0200 
and  0230  to  0330  MST.   Their  height  levels  nearly  coincide  with  the  stratifi- 
cations in  temperature  (Fig.  3.4).   After  0330  MST,  the  low-level  jet  cen- 
tered at  mid- tower  level  expanded  to  a  broad  maximum  as  values  of  wind  speed 
near  the  surface  diminished  gradually.   After  0500  MST,  there  were  two  wind 
speed  maximas  at  100  and  200  m,  at  nearly  the  same  levels  as  the  inversions. 
These  observations  are  in  agreement  with  the  statistical  results  of  Li  et  al, 
(1982b),  which  indicate  a  strong  relationship  between  multiple  wind  speed 
maxima  and  the  existence  of  multiple  inversion  layers.   The  heights  of  the 
wind  speed  maxima  roughly  coincide  with  the  tops  of  the  inversions. 

Figure  3.6  shows  the  time-height  cross  section  of  specific  humidity  for 
the  same  time  as  the  temperature  field  in  Fig.  3.5.   It  is  worth  noting  that 
the  distribution  of  the  isolines  of  specific  humidity  is  similar  to  the 
distribution  of  isotherms  (Fig.  3.4).   After  0230  MST,  the  moisture  field 
also  developed  multiple  layers.   However,  the  height  level  of  the  maximum 
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Figure  3.6. — The  time-height  cross  section  of  humidity  on  24  April  1981, 
Solid  lines  are  isolines  of  specific  humidity  at  0.1-g/kg  intervals; 
thick  dashed  lines  indicate  maximum  values  of  specific  humidity. 


value  of  moisture  is  below  the  top  of  the  inversion.   The  locations  of  the 
strong  moisture  gradients,  on  the  other  hand,  very  nearly  correspond  to  those 
of  the  temperature  gradients.   Radiative-transfer  relations  show  that  the 
greater  the  moisture  content  in  the  inversion  layer  near  the  surface,  the 
stronger  the  radiative  cooling  at  the  top  of  the  inversion.   Cooling  rates  in 
the  presence  of  moisture  will  be  discussed  in  Sec.  3.5. 

3.4   EDDY  KINETIC  ENERGY  AND  THE  DISTRIBUTION  OF  TURBULENT  FLUXES  OF 
HEAT  AND  MOMENTUM 

The  vertical  distributions  of  Richardson  number  (Ri),  horizontal  eddy 

kinetic  energy  (e,),  variances  (w  and  6  ),  and  fluxes  (wu,  wv,  and  wG)  for 
the  10-min  period  from  0250  to  0300  MST  are  calculated  using  the  observa- 
tional data  at  eight  fixed  levels  on  the  BAO  tower.   [Overbars  indicate  time 
average  (10  min);  u,  v,  and  w  are  the  longitudinal,  lateral,  and  vertical 
fluctuations  of  the  wind;  and  9  is  the  temperature  fluctuation.] 
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The  turbulent  fluxes  and  energy  were  calculated  by  integrating  spectral 
and  cospectral  estimates  in  the  frequency  band  0.02-5  Hz.   The  results  are 
presented  in  Fig.  3.7  and  in  Table  3.1.   Ri  and  e,  were  calculated  using  the 
following  relations: 


eh  =  1/2 (u2  +  v2)  , 


and 


Ri  = 


=  £_ 


v 


where  g  is  the  acceleration  of  gravity,  6  is  the  mean  virtual  temperature, 
U  and  V  are  the  mean  horizontal  wind  components  along  the  E-W  and  N-S  axes, 
and  A  denotes  an  increment  over  a  finite  height  interval  Az. 

The  Richardson  number  achieves  its  maximum  value  at  approximately  the 
jet  level  and  has  another  maximum  near  the  center  of  the  inversion  (Fig.  3.7) 
Therefore,  the  Richardson  number  also  has  a  two-layer  structure.   Its  struc- 
ture is  in  agreement  with  the  acoustic  record  of  Fig.  3.2.   The  observations 
of  Mahrt  et  al.  (1979),  based  on  aircraft  soundings  at  Buckley  Field,  Colo., 
in  September  1975,  also  show  that  Richardson  numbers  reach  a  maximum  and 
turbulence  levels  a  minimum  near  the  core  of  a  jet.   The  profiles  of  e,  and 
Ri  in  Fig.  3.7  show  the  same  behavior.   Between  50  and  150  m  the  values  of  Ri 


v'V 


I    I    I    I    I    I    I    I    I 

WV  W0 


V 


8  12 


02(1O"3  °C2) 
ehandw?(10"3m2/s2) 


I  I  I    /\^.    I     I      I     I     I     I      I  T     I  yT' 


I   I  I   I  I   I  I  I  I 


16  -8 


8  16  24 


w0(1O       °Cm/s) 
wuandwv(10~    m2/s2) 


Figure   3.7. — The   distributions   of  Ri,    6    ,   w   ,    e,  ,   wv,   w6 ,    and  wu  with 
height   at  0250  MST  on  24  April   1981. 
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are  close  to  its  subcritical  value  (Ri  =  0.25);  between  250  and  300  m, 
Ri  <  0.25.   These  regions  have  large  e,  values.   That  the  temperature  field 
is  uniform  between  the  two  inversion  layers,  combined  with  the  low  Ri  and 
high  e  values  between  the  layers,  indicates  that  this  region  is  well  mixed 

2  - 

in  the  vertical.   The  distribution  of  w  is  similar  to  that  of  e,  .   Andre 

h 

et  al.  (1978)  gave  the  following  interpretation  for  the  vertical  structure  of 

nocturnal  turbulence:  the  horizontal  components  u  and  v  are  generated  by 

wind  shear,  3u/8z  and  8v/8z;  their  transformation  into  the  vertical  component 

~~2 

w  by  pressure  fluctuations  would  tend  to  give  rise  to  a  strong  negative  heat 

flux. 


The  results  in  Fig.  3.7  indicate  that  wu  is  maximum  and  positive  near 
the  jet  core.   Therefore,  the  vertical  flux  of  longitudinal  momentum  is 
upward.   This  profile  is  consistent  with  the  upward  movement  of  the  jet  core 
between  0250  and  0300  MST  indicated  in  the  wind  speed  cross  section  in 
Fig.  3.5.   However,  wv  is  negative  above  the  jet  core  and  positive  below,  so 
that  the  vertical  transport  of  the  lateral  momentum  is  into  the  jet.   Fig- 
ure 3.7  shows  that  there  is  a  general  tendency  for  the  temperature  variance 

6  to  decrease  slightly  with  height  above  about  50  m. 


3.5   TEMPERATURE  BUDGET  PROFILE 

To  help  understand  the  maintenance  or  breakup  of  multiple  layers  in  the 
stable  nocturnal  atmosphere,  it  is  important  to  study  the  different  factors 
contributing  to  variations  in  the  temperature  profile.   The  factors  include 
radiative  flux  divergence,  turbulence  heat  flux  divergence,  advection,  and 
subsidence. 

The  profile  of  nocturnal  radiative  cooling  with  height  is  complex  when 
temperature  and  humidity  vary  with  height.   To  evaluate  this,  the  Brooks 
(1950)  tabular  method  is  used  to  calculate  nocturnal  radiative  cooling  rates 
with  the  BAO  profile  data.   The  derived  formula  is  the  following: 
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where  q  is  the  specific  humidity;  K  the  pressure  correction  factor;  C   the 
specific  heat  of  air  at  constant  pressure;  co  the  optical  thickness  of  the 
radiating  substance;  a  the  Stefan's  constant;  T  the  absolute  temperature;  £ 
the  emissivity;  co  the  upper  boundary  of  the  radiating  substance;  and  the 
arrows  indicate  the  direction  from  which  radiation  reaches  the  reference 
level.   Subscript  E  denotes  the  value  at  the  top  of  the  radiating  layer. 

4 
If  the  integrals  are  simplified  by  considering  3(CJT.  )/3co  constant  within 

4 
each  integration  interval  and  written  approximately  as  A(aT   )/Aco,  (3.1)  can  be 

written  as 


II  BE 
3t  "  C 
P 


[l/^Aft),^-  f/<>.(Il].  (3-2) 

Ln=l         \   /  m=l         \   /J 


where  3S  /3co  =  OT   (3e  /3co)  is  computed  at  700  mb.   The  summation  over  n  is 
J!        At, 

for  layers  above  the  reference  level  into  which  the  atmosphere  has  been  di- 
vided, and  the  summation  over  m  is  for  layers  below  the  reference  level. 

4 
A(aT   )  is  zero  at  the  ground  surface. 


The  emissivity  curve  3e/3co  was  taken  from  Brooks  (1950,  Table  C) ,  which 
ignores  CO  absorption  (assumed  negligible) .   Brooks  calculates  the  water 
vapor  absorption  from  laboratory  measurements  using  isothermal  radiation  and 
small  path  lengths.   The  vertical  integration  of  (3.2)  begins  at  the  2-m 
level. 

The  variation  of  radiative  cooling  rate  with  height  has  been  calculated 
using  (3.2).   Figure  3.8  shows  the  results  computed  from  the  carriage  profile 
that  began  at  0250  MST.   Because  the  surface  inversion  below  22  m  is  very 
strong  and  moist,  the  radiative  cooling  rate  is  very  strong  (-0.47°C/10  min) 
at  the  top  of  this  layer.   This  distribution  of  temperature  and  humidity  also 
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Figure  3.8. — Profiles  of  specific  humidity  (q),  temperature  (T) ,  total 

cooling  ( ),  radiative  cooling  ( — — ) ,  residual  cooling  (•••"),  and 

turbulent  cooling  ( )  for  the  case  on  24  April  1981. 


accounts  for  the  strong  warming  in  the  layer  nearest  the  surface.   In  addi- 
tion, the  cooling  rate  is  larger  at  the  levels  of  the  temperature  maxima. 
The  cooling  rate  maximum  at  the  top  of  the  surface  inversion  will  cause  the 
inversion  strength  below  to  weaken  with  time,  and  the  temperature  gradient 
above  to  strengthen.   The  net  effect  of  such  cooling  may  be  to  maintain  any 
elevated  inversion  present  (Staley,  1965). 

We  can  study  the  relative  importance  of  radiative  cooling  in  the  tempera- 
ture budget  profile  for  this  case  using  the  relation 


AT  A  — 

j—   =  (Radiative  cooling)  -  -r 1-  (Residual  cooling) 


(3.3) 


AT  /At  is  the  total  temperature  change  during  the  10-min  interval;  the  radia- 
tive cooling  is  calculated  from  (3.2);  the  turbulent  cooling  rates  -Aw6/Az 
are  calculated  from  the  results  in  Table  3.1;  and  the  residual  term  contains 
the  horizontal  advection  of  temperature  and  subsidence.  The  divergence  of  the 
horizontal  heat  flux  is  presumed  to  be  small. 

Table  3.2  presents  the  results  for  height  intervals  between  successive 
tower  levels,  and  Fig.  3.8  the  plots  of  profiles  of  all  the  terms  in  the  heat 
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Table  3.2. — Temperature  budget  (°C/10  min)  for  the 
case  of  24  April  1981 


z 
(m) 


AT 
Ri  At 

(interpolated) 


Turbulent 
cooling 

-Aw6/Az 

(x  10"1) 


Radiative 
cooling 


Residual 
cooling 


16 

0.98 

-1.10 

-0.50 

-0.17 

-0.88 

36 

2.24 

-0.28 

+0.08 

-0.35 

+0.06 

75 

0.38 

+0.70 

+0.14 

-0.08 

+0.77 

125 

0.49 

+0.02 

0 

-0.06 

+0.08 

175 

4.37 

0 

-0.02 

0 

0 

225 

0.36 

-1.43 

-0.10 

-0.09 

-1.33 

275 

0.15 

-0.60 

+0.04 

-0.05 

-0.55 

budget.   We  see  that  the  terms  for  radiative  and  turbulent  cooling  are  much 
smaller  than  the  term  for  residual  cooling  at  all  levels.   The  horizontal 
advection  dominates  the  residual  term,  which  is  very  likely  a  result  of  the 
rapid  change  in  the  height  of  the  inversion  seen  in  Fig.  3.4. 

The  relatively  small  turbulent  cooling  rates  presented  in  Table  3.2  agree 
with  the  model  results  of  Andre  et  al.  (1978)  and  Yamada  and  Mellor  (1975). 
Other  researchers  have  calculated  layer-averaged  temperature  budgets  in  the 
nocturnal  boundary  layer  for  single-layer  cases  (Garratt  and  Brost,  1981; 
Andre  and  Mahrt,  1982).   However,  their  calculations  were  performed  in  the 
nocturnal  boundary  layer,  which  is  only  a  part  of  the  nocturnal  lower  atmo- 
sphere.  Within  the  boundary  layer  the  temperature  fluxes  were  more  than  an 
order  of  magnitude  larger  than  the  case  presented  here.   Therefore,  it  is 
difficult  to  compare  their  results  with  those  in  Table  3.2. 

The  importance  of  advection  and/or  horizontal  divergence  of  heat  flux  is 
clearly  shown  by  the  large  values  of  the  residual  cooling  listed  in  Table  3.2. 
The  two-dimensional  models  cited  above  certainly  do  not  apply  to  this  fairly 
typical  case  studied  at  the  BAO,  which  is  located  in  rolling  terrain.   Rea- 
listic predictions  for  stable  boundary  layers  must  come  from  detailed  three- 
dimensional  models. 
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3.6   CONCLUSIONS 

The  above  analysis  indicates  that  the  development  and  the  structure  of 
the  stable  nocturnal  lower  atmosphere  are  fairly  complicated.   Multiple 
layers  are  observed  over  both  uniform  and  complex  terrain  and  are  most  readily 
apparent  in  acoustic  sounder  records.   They  are  also  apparent  in  plots  of  the 
wind  field,  the  temperature  field,  and  the  moisture  field.   Atmospheric  and 
topographical  factors  contributing  to  such  layering  are  many,  but  horizontal 
advection  is  involved  in  some  way  in  their  formation. 

Because  the  Richardson  number  achieves  its  maximum  value  near  the  center 
of  a  low- level  jet,  there  appears  to  be  a  rough  correspondence  between  the 
profiles  of  Richardson  number  and  mean  wind  speed.   In  the  region  between  its 
maxima,  Ri  tends  to  be  subcritical  and  e,  to  be  maximum.   In  the  same  region, 
the  temperature  field  is  fairly  uniform,  suggesting  the  region  between  two 
elevated  layers  is  well  mixed.   The  turbulent  vertical  fluxes  of  heat  within 
elevated  layers  and  near  the  ground  surface  are  negative,  but  the  vertical 
distribution  of  momentum  flux  is  very  complex. 

The  vertical  distribution  of  radiative  cooling  rate  in  a  nocturnal  lower 
atmosphere  shows  warming  nearest  the  surface  and  cooling  aloft.   Cooling  was 
most  pronounced  at  the  top  of  the  strong  inversion  layer  near  the  surface 
(z  *  22  m).   The  total  heat  budget  for  the  multiple-layer  case  showed  that  the 
turbulent  cooling  rates  were  negligible  compared  with  the  radiative  cooling, 
and  both  were  much  smaller  than  horizontal  advection  and/or  subsidence. 
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ABSTRACT.   The  phase  speed  and  azimuth  of  gravity  waves  in  the 
planetary  boundary  layer  are  calculated  using  an  optical  tri- 
angle.  The  results  are  compared  with  those  calculated  using  a 
microbarograph  array.   The  two  methods  are  in  essential  agree- 
ment. 


4 . 1   INTRODUCTION 

Gravity  waves,  which  are  often  synoptically  generated,  are  frequently 
observed  in  the  nocturnal  stable  boundary  layer.   Because  of  the  rapid  develop- 
ment of  atmospheric  remote  sensing  techniques  in  recent  years,  the  temporal 
and  spatial  structure  of  gravity  waves  can  be  displayed  and  studied.   Gossard 
and  Richter  (1970),  Emmanuel  et  al.  (1972),  and  Hooke  et  al.  (1973),  among 
others,  have  used  sodars  and  FM-CW  radars  to  monitor  the  physical  characteris- 
tics of  gravity  waves  and  to  study  the  interaction  between  the  waves  and 
turbulence.   These  authors  have  also  studied  the  relationship  between  internal 
gravity  waves  and  synoptic  and  mesoscale  phenomena.   The  studies  have  con- 
sidered the  causes  and  the  effects  of  the  waves  on  mesoscale  activity. 
The  phase  speed  and  the  azimuth  of  gravity  waves  are  very  important  to  help 
understand  their  source.   In  these  studies,  the  phase  speed  and  the  azimuth  of 
gravity  waves  generally  were  measured  using  the  time  series  data  of  three  or 
more  microbarographs.   In  addition,  Gossard  and  Munk  (1954)  used  single-point 
data  for  pressure  and  wind  speed  to  compute  the  phase  speed.   Recently,  Eymard 
and  Weill  (1979)  have  measured  the  phase  speed  and  the  azimuth  of  gravity 
waves  (along  with  other  characteristics)  using  a  three-antenna  Doppler  sodar. 
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The  following  describes  the  calculations  of  some  basic  characteristics  of 
the  gravity  waves  utilizing  an  optical  triangular  array  surrounding  the  300-m 
meteorological  tower  at  the  Boulder  Atmospheric  Observatory  (BAO)  in  Colorado. 
The  wind  speed  measured  by  the  optical  systems  is  an  average  along  the  path 
(Lawrence  et  al. ,  1972).   Therefore,  high  frequencies  are  automatically  fil- 
tered.  Because  gravity  waves  are  low-frequency  phenomena  by  atmospheric 
boundary  layer  standards,  the  low  pass  filtering  of  the  optical  sensors  may 
improve  phase  speed  and  direction  estimations  when  compared  with  an  array  of 
microbarograph  sensors.   Measuring  the  horizontal  divergence  of  the  wind  speed 
with  an  optical  triangle  (Kjelaas  and  Ochs,  1974;  Tsay  et  al.,  1980)  has 
proved  effective.   With  this  information  in  mind,  measuring  the  phase  speed 
and  azimuth  of  gravity  waves  with  optical  methods  appears  to  be  a  reasonable 
approach. 


4.2  MEASUREMENTS  AND  COMPUTATIONS 

Three  optical  systems  have  been  operating  at  BAO.   The  optical  paths  form 
an  equilateral  triangle,  450  m  on  a  side  (Fig.  4.1).   Because  the  optical 
sensors  measure  the  mean  wind  perpendicular  to  each  path  and  positive  outward 
from  the  triangle  center  (Lawrence  et  al. ,  1972),  the  sum  of  the  readings  of 
all  three  paths  represents  the  average  horizontal  wind  speed  out  from  the 
triangle.   The  individual  wind  readings  for  each  optical  path,  smoothed  with  a 
10-s  running  average  along  with  the  scalar  sum  of  all  three,  were  recorded  on 
digital  tape.   For  the  three  optical  systems,  the  horizontal  area-average  wind 
speed  and  the  wind  direction  can  be  expressed  as 


V  =  (VE  +  VN  j 


(4.1) 


and  v 

G  =  tan  -1  rp  +  180°  +  4°,  (4.2) 

N 

where  V  and  V  are  the  wind  components  oriented  to  the  east  and  north  respec- 

tively,  VE  =  A-  (VSE  -  VSH),  and  ?„  =  §  T„  -  \   (VgE  +  ?„).   VN>  V^,  and  V,,,, 

are  the  wind  components,  transverse  to  the  north,  southeast,  and  southwest 
triangle  sides  repectively.   In  Eq.  (4.2),  adding  180°  makes  the  northerly  winds 
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Figure  4.1. — Schematic  of  the  optical  triangle,  the  three  microbarographs 
used  in  this  analysis  (x) ,  and  the  BAO  tower  (•)•   V.  =  wind  component 

transverse  to  the  triangle  sides;  V!  =  actual  wind  speed  for  each  tri- 
angle side;  a  =  angle  between  V.  and  V!. 


represent  a  0°  direction,  and  the  4°  is  an  adjustment  for  the  difference 
between  the  north  path  of  the  optical  triangle  and  the  true  east-west  direc- 
tion.  The  wind  directions  9  correspond  to  the  standard  meteorological  defi- 
nition. 

The  theoretical  wind-weighting  function  for  a  circular  transmitter  and 
receiver  with  circular  apertures,  3.0  and  1.95  Fresnel  zones  in  diameter, 
respectively,  has  its  peak  at  the  center  of  the  path  as  shown  in  Fig.  4.2  (Ochs 
et  al. ,  1976).   The  curve  of  the  wind-weighting  function  also  shows  the  mini- 
mum value  near  the  transmitter  and  receiver  along  the  optical  path.   There- 
fore, the  average  wind  speed  measured  by  optical  systems  is  strongly  weighted 
toward  the  center  of  the  path,  and  it  is  assumed  that  a  new  triangle  can  be 
constructed  by  joining  the  center  points  of  each  optical  path. 

For  each  optical  path,  the  actual  wind  speed  at  the  center  point  of  the 
path  (Fig.  4.1)  can  be  written 
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Figure  4.2. — Experimental  wind-weighting  function  obtained  by  comparison 
of  the  optical  measurement  over  a  500-m  path  with  five  pairs  of  anemom- 
eters arranged  to  measure  the  horizontal  cross-wind  component  (from 
Ochs  et  al. ,  1976).   T  =  transmitter;  R  =  receiver. 


V!  =  V./cosa.  , 


(4.3) 


where  V.  is  the  wind  component  transverse  to  the  triangle  sides,  V'.  is  the 
actual  wind  speed  for  each  triangle  side,  the  subscript  i  (=  1,2,3)  represents 
the  index  of  each  side  of  the  optical  triangle,  and  a.  is  the  angle  between 
the  average  wind  direction  and  the  wind  component  transverse  to  each  triangle 
side.   Because  the  average  wind  speed  and  wind  direction  change  with  time,  the 
angle  a.  is  also  time  variable.   The  angles  made  with  each  side  of  the  optical 
triangle  may  be  instantaneously  different.   Therefore,  with  gravity  wave 
propagation,  the  phase  speed  and  direction  appear  as  a  variation  in  the  time 
series  of  the  actual  wind  speed  for  the  three  sides.   The  propagation  of  a 
wave  across  the  three  separate  optical  systems  can  be  distinguished  by  the 
time  delays  between  the  three  time  series  of  wind  speeds  for  pairs  of  optical 
stations.   We  can  then  obtain  the  phase  difference  between  the  pairs  of 
optical  stations  for  each  frequency. 


If  the  phase  speed  and  the  propagation  direction  of  gravity  waves  are 
represented  as  c  and  y   respectively,  then 


96 


Wave  train 


/ 
/ 
/ 


Direction  of 
gravity  wave 


Figure  4.3. — Schematic  of  the  geometry  to  determine  the  time  delay  with 

the  three  separate  optical  systems,  y  =  propagation  direction  of  gravity 

waves;  <J)  =  azimuth;  £  =  length  between  the  central  points  of  the  sides 
of  the  optical  triangle. 


and 


cosy  =  c  t12/£ 


cos(60°  -  y)  =  c  t13/£, 


(4.4) 
(4.5) 


where  i   is  the  length  between  the  central  points  of  the  sides  of  the  optical 
triangle  and  t  _  and  t  _  are  the  time  delays  from  vertex  1  to  vertex  2  and  to 
vertex  3  respectively  (Fig.  4.3).   The  corners  of  the  triangle  in  Fig.  4.3  are 
at  the  center  points  of  each  optical  path. 


From  Eqs.  (4.4)  and  (4.5),  c  and  Y  can  then  be  written 


c  =  (£/t12)  VI  +  tan"  Y 


A 


Y  =  tan 


'fefel 


(4.6) 
(4.7) 
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with 


where  <j)  is  the  azimuth. 


<J>  =  210°  -  y  , 


(4.8) 


To  calculate  c  and  (J),  we  must  derive  the  time  delays  t  „  and  t  „.   If 
the  time  series  of  the  wind  speed  is  V (t)  for  vertex  1  and  V'(t)  for  vertex 
2,  the  power  spectrum  is  obtained  by  Fourier-transforming  the  auto-covariance 
function,  and  the  cross  spectrum  is  obtained  by  transforming  the  cross-covari- 
ance  function  between  the  two  time  series  (Gossard  and  Hooke,  1975;  Brigham, 
1974).   Therefore,  the  form  of  the  cross-power  spectrum  for  V  (t),  V'(t)  can 
be  written  as 


oo 

E(U>)  =  /  [/  V'(t)  Vj(t  +  T)  dl 


-1C0T  , 

e     dT   , 


(4.9) 


where  T  is  the  time  delay  between  time  series,  and  go  =  2iTf,  with  f  the  fre- 
quency. 


If  a  =  t  +  T,  (4.9)  becomes 

oo 

E(a))  =  /  V|(t)  eiaJt  H(a))  dt  , 


(4.10) 


where 


H(a>)  =  /  Vl(a) 


-icoa  , 
e     da 


(4.11) 


Finally,  from  Eq.  (4.10),  the  form  of  the  cross-power  spectrum  can  be 
written  as 


E(oj)  =  H(cj)  [R(w)  +  il(0))]  , 


(4.12) 


where 


and 


R(co)  =  /  V'(t)  cos(oJt)  dt 


I(oo)  -  /  V'(t)  sin(wt)  dt 


(4.13) 


(4.14) 
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The  phase  relationship  between  frequencies  for  two  cross-spectral  analyzed 
time  series  is  given  by 

tan  ojt  =  I(u))/R(oo)  ,  (4.15) 

where  T  =  t   ,  in  this  case.   The  time  delay  t..  _  can  be  similarly  calculated. 
The  phase  speed  and  direction  of  the  gravity  wave  can  be  obtained  by  Fourier- 
transforming  the  auto-covariance  function  (from  the  time  series  of  the  wind 
speed)  using  Eqs.  (4.12),  (4.13),  and  (4.14)  and  solving  Eqs.  (4.6),  (4.8), 
and  (4.15). 


4.3  RESULTS 

Four  gravity  wave  cases  are  presented.   The  first,  between  0000  and 
0700  MST  on  13  November  1981  was  embedded  in  strong  stable  stratification  at 
the  lower  levels.   The  gravity  wave  train  is  shown  in  the  acoustic  record  pre- 
sented in  Fig.  4.4.   The  phase  speed  and  azimuth  of  the  gravity  waves  have 
been  calculated  from  0210  to  0330  and  from  0530  to  0650  MST.   We  first  discuss 
the  gravity  wave  occurring  between  0210  and  0330  MST.   The  optical  wind  sensor 
was  operating  with  a  f ive-microbarograph  array  near  the  BAO  tower.   The  three 
microbarographs  closest  to  the  tower  were  used  for  the  analysis.   Their  loca- 
tions relative  to  the  optical  triangle  are  presented  in  Fig.  4.1.   The  time 
series  of  the  three  microbarographs  are  shown  in  Fig.  4.5a.   The  coherence  of 
the  traces  is  quite  high.   The  time  series  of  the  three  optical  systems  are 
shown  in  Fig.  4.5b.   Similarly,  their  coherence  is  very  good.   The  mean  spec- 
trum of  the  three  optical  wind  measurements  is  presented  in  Fig.  4.6a.   The 
major  peak  is  at  a  34.1-min  period.   Two  minor  peaks  are  at  19-  and  14.2-min 
periods.   The  phase  speed  and  azimuth  are  14.4  m/s  and  200°  respectively.   The 
mean  spectrum  is  also  calculated  for  the  same  time  period  (Fig.  4.6b)  using 
the  three  microbarographs.   The  major  peak  appears  at  a  28.4-min  period  with 
minor  peaks  at  19-  and  14.2-min  periods.   The  phase  speed  and  azimuth  of  the 
major  peak  are  17.7  m/s  and  225°  respectively. 

For  the  0530  to  0650  MST  time  period  (Fig.  4.4),  the  two  sets  of  three 
time  series  (not  shown)  are  also  very  coherent.   Figures  4.7a,b  indicate  a 
similarity  between  the  optical  wind  and  microbarograph  mean  spectra  with 
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Figure  4.5a. — The  time  series  of  the  three  microbarographs  between  0210 
and  0335  MST  on  13  November  1981.  The  vertical  scale  between  two  lines 
is  80  yb/s,  and  the  curves  are  relative  to  a  zero  mean. 
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Figure  4.5b. — The  time  series  for  three  separate  optical  systems  from 
0210  to  0335  MST  on  13  November  1981.   The  vertical  scale  between  two 
lines  is  3  m/s,  and  the  numbers  on  the  far  left  are  the  individual  means. 
V  ,  V   ,  and  VOIT  are  the  wind  components,  transverse  to  the  north,  south- 

east,  and  southwest  triangle  sides,  respectively. 
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Figure  4.6a. — The  mean  spectrum  calculated  from  the  optical  triangle 
data  from  0210  to  0330  MST  on  13  November  1981. 
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Figure  4.6b. — The  mean  spectrum  calculated  from  the  three  microbarographs 
from  0210  to  0330  MST  on  13  November  1981. 
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Figure  4.7a. — The  mean  spectrum  calculated  from  the  optical  triangle  data 
from  0530  to  0650  MST  on  13  November  1981. 
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Figure  4.7b. — The  mean  spectrum  calculated  from  the  three  micorbarographs 
from  0530  to  0650  MST  on  13  November  1981. 
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Table  4.1. — Periods  of  spectral  peaks,  phase  speeds,  and  azimuths 
calculated  from  the  optical  triangle  (0)  and  microbarographs  (M) 

for  the  four  cases  discussed 


Bandwidth 

Periods  of 

of  wave 

Phase 

Date  and  time 

Array 

peaks 

analysis 

speed 

Azimuth 

(MST) 

(min) 

(min) 

(m/s) 

(deg) 

13  Nov.  1981 

0 

34.1, 

19,  14.2 

3.4-85.3 

14.4 

200 

0210-0330 

M 

28.4, 

19,  14.2 

3.4-85.3 

17.7 

225 

13  Nov.  1981 

0 

28.4, 

21.3 

3.4-34.1 

10.0 

254 

0530-0650 

M 

24.4, 

17.1 

3.4-34.1 

10.0 

270 

27  Oct.  1981 

0 

34.1, 

24.4,  17. 

,1 

3.4-85.3 

12.0 

360 

1940-2100 

M 

34.1, 

24.4,  13. 

,1 

3.4-85.3 

10.0 

360 

28  Oct.  1981 

0 

21.3, 

14.2 

3.4-85.3 

11.6 

296 

0530-0650 

M 

21.3, 

14.2,  12. 

,2 

3.4-85.3 

10.0 

270 

respect  to  the  major  peaks.   The  phase  speed  from  both  the  optical  array  and 
microbarograph  is  10  m/s,  and  the  azimuths  are  254°  and  270°  respectively. 

The  two  other  gravity  wave  events  analyzed  occurred  on  27  and  28  October 
1981.   A  summary  of  results  for  these  periods  is  given  in  Table  4.1  along 
with  results  for  events  of  13  November  described  above.   The  table  indicates 
that  the  agreement  of  the  phase  speeds  and  directions  between  the  two  methods 
is  quite  good. 

One  disadvantage  in  using  the  optical  triangle  for  gravity  wave  studies 
can  be  seen  from  Eq.  (4.3)  and  Fig.  4.1.  When  the  actual  wind  V.  approaches  a 
parallel  direction  with  one  of  the  optical  paths  (a.  •*■   90°),  the  wind  becomes 
indeterminate  for  that  path.   In  this  situation,  a  data  spike  will  appear  in 
the  V!  data  stream  for  that  path.   The  spikes  can  be  removed  by  interpolation 
using  the  acceptable  points  on  either  side. 

As  mentioned  earlier,  the  spatial  filtering  of  the  optical  array  can  be 
an  advantage.  One  advantage  to  this  built-in  filtering  is  in  the  use  of  the 
impedance  relation  (Gossard  and  Munk,  1954), 


HI   =   P'(z> u   (z) 

k   p  (z)u' (z)    okVZ;  » 


(4.16) 
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where  both  the  background  wind  U   and  the  wave-associated  departure  of  the 
wind  u1  can  be  measured  with  the  optical  triangle.   In  Eq.  (4.16),  p'  is  the 
wave-perturbed  pressure,  p  the  mean  density,  and  k  the  wave  number.   Often, 
using  a  time  series  of  the  wind  from  a  single  point  sensor  presents  a  diffi- 
cult filtering  problem  in  ascertaining  a  proper  u' .   The  optical  triangle 
winds  alleviate  many  of  the  filtering  problems  (E.E.  Gossard;  NOAA/ERL/WPL; 
private  communication). 

However,  the  spatial  filtering  of  the  optical  winds  can  also  present  a 
limitation  not  associated  with  the  microbarograph  technique.   Waves  with  wave 
lengths  equal  to  or  smaller  than  the  averaging  distance  (450  m)  will  be  fil- 
tered out  completely  using  the  optical  technique.   Wave  lengths  slightly 
longer  than  450  m  will  be  filtered  considerably  depending  on  the  filter  func- 
tion of  the  unit.   For  the  major  peaks  analyzed  in  this  study,  all  of  the 
wavelengths  are  significantly  longer  than  450  m.   A  comparison  of  the  optical 
wind  and  microbarograph  spectra  in  Figs.  4.6  and  4.7  shows  how  the  spatial 
filtering  with  the  optical  triangle  affects  the  higher  frequencies. 


4.4   CONCLUSIONS 

The  phase  speeds  and  directions  of  gravity  waves  are  possible  to  calcu- 
late, within  the  limits  of  the  spatial  averaging  and  nonparallel  wind  direc- 
tions, using  an  optical  triangle  wind-measuring  system.   The  results  agree 
quite  well  with  those  from  a  microbarograph  array.   Some  of  the  differences 
between  the  optical  wind  spectra  and  microbarograph  spectra  can  be  explained 
by  optical  wind  spatial  filtering.   The  advantages  of  the  optical  triangle 
over  a  microbarograph  array  include  the  inherent  spatial  averaging  of  the 
optical  triangle.   The  spatial  filtering  can  be  a  disadvantage  for  wavelengths 
near  the  pathlength. 

A  discussion  of  the  connection  between  wave-induced  pressure  fluctuations 
and  wave-induced  wind  fluctuations  near  the  ground  under  statically  stable 
conditions  is  beyond  the  scope  of  this  work.   However,  consideration  of  such  a 
connection  is  important  to  assess  the  ultimate  applicability  of  optical  wind 
sensing  to  atmospheric  gravity  wave  measurement. 
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5.   RICHARDSON  NUMBER  COMPUTATIONS  IN  THE  PLANETARY  BOUNDARY  LAYER 

R.J.  Zamora* 
Department  of  Meteorology 
Metropolitan  State  College 
Denver,  Colorado  80204 

ABSTRACT.   This  paper  focuses  on  a  method  of  calculating  the 
gradient  Richardson  number  Ri  using  data  gathered  at  the  Boulder 
Atmospheric  Observatory,  a  300-m  meteorological  tower  located 
30  km  north  of  Denver,  Colorado,  and  presents  a  case  study  of 
Richardson  number  (Ri)  behavior  at  the  75-m  level  of  the  tower 
during  an  active  acoustic-gravity  wave  event.   Ri  is  found  to  be 
modulated  by  the  gravity  wave  activity  in  addition  to  longer 
period  oscillations  that  might  be  a  result  of  inertial  wave 
activity. 


5.1   INTRODUCTION 

Recent  experiments  conducted  by  NOAA/ERL/WPL  have  focused  on  the  struc- 
ture of  the  stable  planetary  boundary  layer  (PBL) .   One  of  the  more  important 
parameters  used  to  quantitatively  estimate  the  dynamic  stability  of  such  a 
fluid  environment  is  the  gradient  Richardson  number  Ri.   This  has  long  been 
regarded  as  an  important  criterion  for  the  onset  of  turbulence  in  a  stably 
stratified  fluid.   Fluid  parcels  can  be  moved  vertically  against  the  force  of 
gravity  if  sufficient  kinetic  energy  is  available  in  the  velocity  field. 

Ri  represents  the  ratio  of  work  needed  to  bring  about  the  exchange  to  the 
kinetic  energy  present  to  do  the  work.   In  general  if  Ri  >  0.25  through  the 
fluid,  the  fluid  is  assured  to  be  stable  everywhere.   These  results  were 
obtained  by  Miles  (1961)  and  Howard  (1961)  for  incompressible  fluids,  and  the 
proof  was  later  extended  to  compressible  fluids  by  Chimonas  (1970). 

Although  the  Richardson  number  is  of  great  theoretical  importance,  it  is 
a  difficult  parameter  to  calculate.  In  its  usual  form  for  atmospheric  appli- 
cations (Gossard  and  Hooke,  1975), 


*  Present  affiliation:  NOAA/ERL/Wave  Propagation  Laboratory,  Boulder, 
Colorado  80303. 
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with 


N2 
Ri  =  -^—     *  (5.D 

dz 


^-lfz        ,  (5.2) 


where  g  is  the  gravitational  acceleration,  N  the  Brunt-Vaisala  frequency,  6 
the  potential  temperature,  z  the  vertical  Cartesian  coordinate,  U  the  hori- 
zontal wind  vector,  and  9  the  mean  potential  temperature  over  the  gradient. 

To  fully  appreciate  the  difficulty  in  calculating  Ri ,  assume  Ri  =  0.25,  a 

•¥■  2 

typical  value  in  a  stable  layer  with  shear,  and  (9U/9z)   =  0.004,  a  moderate 

value.   Then 

|£  =  1.5  x  10"3  °C/m  . 
dz 

These  correspond  to  differences  of  approximately  0.05°C  and  0.2  m/s  over  50-m 
Az.   Although  it  is  possible  to  measure  these  differences  with  present  BAO 
instrumentation,  it  must  be  noted  that  absolute  calibration  errors  can  bias 
the  calculations,  and  aliasing  due  to  the  finite  difference  scheme  can  yield 
poor  results.   Thus,  in  cases  where  the  wind  shears  and  temperature  gradients 
are  weak,  the  limitations  of  instrumentation  are  most  likely  to  be  apparent  in 
the  calculations. 

In  calculating  Ri  it  is  commonly  assumed  that  the  atmosphere  is  dry.   If 
the  gradient  of  0  in  the  layer  over  which  the  Richardson  number  is  calculated 
is  less  than  zero,  the  atmosphere  is  statically  unstable  and  any  perturbation 
in  the  atmosphere  will  grow.   In  this  case  the  Richardson  number  is  of  little 
significance  since  the  energy  contained  in  the  velocity  field  is  not  necessary 
for  vertical  parcel  exchange.   If  the  gradient  of  potential  temperature  is 
greater  than  zero,  a  restoring  force  exists  and  a  displaced  parcel  will  oscil- 
late about  its  equilibrium  position  with  a  frequency  equal  to  N  in  the  absence 
of  viscous  damping.   Since  the  potential  temperature  of  a  parcel  depends  only 
on  pressure  and  temperature,  errors  can  be  introduced  into  Ri  computations  if 
the  density  of  the  parcel  is  assumed  to  depend  only  on  the  dry  air  pressure 
and  temperature  of  the  parcel. 
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In  an  attempt  to  relax  the  assumption  of  a  dry  atmosphere,  the  virtual 
temperature  of  a  parcel  can  be  substituted  into  the  potential  temperature 
equation,  and  in  turn  this  virtual  potential  temperature  can  be  used  to  cal- 
culate a  density-adjusted  static  stability  term  and  thus  adjust  Ri  for  the 
effect  of  moisture  in  the  atmosphere. 

The  preceding  stability  analysis  remains  valid  until  the  atmosphere 
reaches  saturation,  at  which  time  the  parcel  cools  at  the  moist  adiabatic 
lapse  rate.   It  has  been  demonstrated  (Lalas  and  Einaudi,  1973)  that  at  this 
point  a  new  Brunt-Vaisala  frequency  must  be  used  to  describe  parcel  buoyancy 
adequately.   This  paper  does  not  address  the  problem  of  dynamic  stability  in  a 
saturated  atmosphere,  or  after  condensation  has  occurred. 

To  examine  the  behavior  of  the  Richardson  number  in  the  lower  atmosphere 
a  data  set  was  chosen  from  a  period  of  active  acoustic-gravity  wave  propaga- 
tion beginning  on  the  evening  of  6  November  1981  at  about  2000  MST  and  ending 
on  7  November  1981  at  0400  MST  at  the  Boulder  Atmospheric  Observatory  (BAO) . 
BAO  is  located  20  km  east  of  Boulder  on  gently  rolling  terrain. 

5.2   GOVERNING  EQUATIONS 

5.2.1  Dry  Atmosphere 

The  gradient  Richardson  number  is  given  by  (5.1),  where  (5.2)  represents 
the  static  stability  of  the  atmosphere,  and 


8U 


■  ¥f  *  ¥f  • 


(5.3) 


This  expression  for  the  change  in  wind  shear  with  height  takes  into  account 
both  the  directional  and  speed  shears  of  the  wind  vector: 


-*• 


U  =  u£  +  v5   .  (5.4) 

The  parcel  potential  temperature  9  is  given  by 

e  =  i(ig2o\K  (5.5) 


/iooo\K 
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where  K  =  0.286,  T  is  the  Kelvin  temperature  of  the  parcel,  and  p  is  the 
atmospheric  pressure  in  millibars. 

5.2.2  Moist  Atmosphere 


The  density- ad jus ted  Richardson  number  Ri  is 

N2 
Ri     =    ,-zk         ,  (5.6) 


2 


where 


'  ¥) 


i2  =  &-  —i 


N     =  *-    t~       ,  (5.7) 

v       0        dz 
v 


0       +0 
v  V 

e    =  _A_ — l 

v  2 


and  virtual  potential  temperature  0  is  given  by  Fleagle  and  Businger  (1980) 


as 


/iooo\K 

v )  ■ 


ev  =  t*|       -ji 


T  is  given  by  Hess  (1959)  as 


I'-iVrf      •  (5.9) 


where  q,  the  mixing  ratio,  is 


q  =  e  — S—  ,  (5.10) 

p  -  G 


with  £  =  0.622. 


The  vapor  pressure  e  of  the  parcel  is  given  by  Murray  (1967)  as 


a(T  -  273.16) 
e  =  6.1078  exp  I  £ r I  ,  (5.11) 


ra(Td  -  273.16)  "1 
L   Td"b     J    ' 
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Figure  5.1. — 8  and  9     vs.    T-T,   for  constant  values  of  p  and  T. 
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where 

a  =  21.8745584  ) 

>  over  ice, 

b  =  7.66       ) 

a  =  17.2693882   ) 

>  over  water, 

b  =  35.86       ) 

T  is  the  parcel  dewpoint. 

The  use  of  6  in  (5.7)  allows  the  computation  of  a  more  realistic  Richard- 
son number  because  the  static  stability  of  the  atmosphere  is  no  longer  assumed 
to  be  a  function  of  p  and  T  only.   Figure  5.1  illustrates  the  differences 

between  9  and  8   for  increasing  T,  at  constant  values  of  p  and  T.   This  ad- 
v  d 

justment  is  negligible  at  low  temperatures  because  e  is  typically  small. 


5.3  COMPUTATIONAL  SCHEME 

5.3.1   Instrumentation  and  Calibration 

The  BAO  facility  includes  a  300-m  meteorological  tower  instrumented  at 
eight  levels  as  shown  in  Fig.  5.2.   Standard  instrumentation  at  each  level 
includes  fast-response  wind  measurements  by  sonic  anemometers,  slow-response 
wind  measurements  using  R.M.  Young  propeller-vane  (propvane)  anemometers,  mean 
temperatures  using  Hewlett-Packard  quartz  thermometers  in  an  aspirated  shield, 
mean  dewpoints  using  EG&G  110  Dewpoint  Hygrometers,  and  fast-response  tempera- 
tures using  platinum  wire  thermometers.   A  full  description  of  BAO  is  given  by 
Kaimal  and  Gaynor  (1983). 

For  the  purpose  of  Richardson  number  calculation  the  propvane  anemometers 
were  chosen  for  velocity  field  measurements,  and  the  mean  temperature  and 
dewpoint  sensors  were  chosen  to  measure  the  thermodynamic  variables.   These 
choices  were  made  for  operational  reasons.   Since  Ri  involves  the  gradients 
30/3z  and  3U/8z,  the  absolute  calibrations  must  be  as  accurate  as  possible. 
The  propvane  anemometers,  and  mean- temperature  and  mean  dewpoint  sensors, 
offered  the  most  accurate  absolute  calibration.   The  mean-temperature  sensors 
are  periodically  placed  in  a  precision  temperature  bath  and  intercompared. 
The  dewpoint  hygrometers  are  calibrated  with  precision  resistors  and  are 
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Figure  5.2. — BAO  tower,  showing  the  location  of  fixed  levels  and  the 
parameters  measured. 
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cleaned  periodically  to  maintain  precision.   The  propvane  anemometers  are 
calibrated  using  a  constant  rpm  motor  for  speed  and  are  oriented  with  respect 
to  the  supporting  boom  for  azimuth. 

The  accuracies  for  the  respective  measurements  are 
+  0.05°C  for  temperature, 
+  0.5°C  for  dewpoint, 
+0.1  m/s  for  speed, 
+1°  for  azimuth. 


5.3.2   The  Algorithm 

To  calculate  the  Richardson  number,  an  interactive  program,  which  resides 

in  the  permanent  library  of  the  BAO  PDP  11/70  computer,  was  written  by  the 

author.   The  source  code  is  written  in  Fortran  IV.   The  program  allows  the 

input  quantities  to  be  averaged  over  variable  increments  of  time  and  allows 

the  option  of  calculating  Ri  or  Ri  . 

v 

Calculation  of  q,  G,  and  6  requires  knowledge  of  p  at  each  tower  level, 
but  p  is  measured  only  at  the  surface.  Thus,  the  hypsometric  equation  (Hess, 
1959)  is  solved  in  the  form 


P2  =  exp 


I""  ^T  (z2  "  zi>l  +  An  Pl   ,  (5.12) 


—■k  *     * 

where  T  =  (T  +  T~)/2,  z~  and  z  are  the  tower  levels  chosen,  and  p..  = 

■k 

surface  pressure.   Because  T  also  requires  knowledge  of  p  at  each  level, 
(5.12)  cannot  be  solved  directly. 

Instead,  the  mean  temperature  T  is  substituted  for  T  in  (5.12)  and  this 
approximate  equation  is  solved  for  p  at  each  tower  level.   This  first  guess  is 
then  used  to  compute  q  and  T  for  the  tower  levels. 


After  this  step  is  completed,  (5.12)  is  then  solved  again  obtaining  a 
corrected  p  for  levels  one  through  eight.   This  p  is  then  used  in  successive 
calculations. 
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The  partial  derivatives  90/9z,  99  /9z,  and  9U/9z  are  approximated  by 
centered  differences  as 


and 


9U 


9z 


90    . 

2 

1 

9z 

Z2   " 

2i   ' 

0 

-   0 

90 

V 

.      V2 

Vl 

9z 

Z2 

"   Zl 

(5.13) 


(U2  "  Ul)2  +   (V2  "  Vl)2 
(z2  -  Z;L)2 


(z2  -  z^y 


(5.14) 


5.4   ERROR  ANALYSIS 

To  test  the  sensitivity  of  the  algorithm  to  experimental  and  truncation 
error,  a  separate  version  of  the  program  was  created  in  which  the  errors  given 
in  Sec.  5.3.1  were  added  to  the  measured  quantities.   The  differences  between 
the  two  computer  runs  are  as  follows : 


ARi  =  0.007 
v 

ARi  =  0.008  . 


Hence,  the  estimates  of  Ri  and  Ri  are  accurate  to  +  0.007  and  +  0.008  re- 

v  —  — 

spectively.   Much  larger  errors  are  possible  if  the  differences  in  wind  veloc- 
ity and  temperature  over  the  layer  in  question  exceed  the  limits  of  the  sensor, 


5.5  ANALYSIS 

The  echo  patterns  from  an  acoustic  sounder  located  near  the  tower 
(Fig.  5.3)  were  used  to  determine  the  onset  of  wave  activity,  and  the  Richard- 
son number  was  computed  between  the  50-  and  100-m  levels  of  the  tower  starting 
at  2000  MST  on  6  November  1981  and  ending  at  0420  on  7  November  1981,  for  an 
averaging  period  of  5  min  using  the  algorithm  of  Sec.  5.3.2.   A  histogram  of 
the  distribution  of  Ri   for  the  period  is  shown  in  Fig.  5.4. 
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Figure  5.3. — BAO  facility  showing  the  location  of  the  acoustic  sounder 
and  microbarograph. 
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Figure  5.4= — Histogram  of  Ri  at  75  m. 


It  was  hoped  that  since  an  acoustic-gravity  wave  represents  the  response 
of  the  atmosphere  to  accelerations  caused  by  deviations  from  geostrophic 
balance,  one  would  observe  a  periodicity  in  Ri  that  would  represent  the 
attempt  of  the  atmosphere  to  gain  geostrophic  balance  by  transfering  momentum 
into  the  lower  regions  of  the  PBL  through  mechanical  turbulence.   Since  the 
gradient  Richardson  number  relates  the  mean  field  gradients  of  temperature  and 
velocity  to  the  dynamic  stability  of  the  atmosphere  and  the  onset  of  mechani- 
cal turbulence,  any  changes  in  the  mean  fields  of  velocity  and  temperature 
should  be  reflected  in  the  turbulence  distribution  (Neff,  1980). 


As  the  wave  propagates  through  the  tower,  parcels  of  differing  0  will 
alternately  move  past  the  sensors,  and  thus  dQ   /dz   should  vary  with  the  fre- 
quency of  the  wave.   Since  the  vertical  momentum  distribution  is  changing, 
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continuity  requires  that  the  horizontal  momentum  distribution  must  also 

-*■    2 

change.   Thus,  not  only  should  36  /9z  vary  with  time,  but  also  (9U/9z) 

must  vary.   Also,  9p/9t,  at  the  ground,  should  reflect  the  wave  propagation. 


The  evening  of  6  November  and  early  morning  of  7  November  are  character- 
ized by  moderate  static  stability  as  shown  by  the  National  Weather  Service 
soundings  taken  at  Stapleton  Airport  near  Denver,  Colo.  (Fig.  5.5).   In  addi- 
tion the  evening  wind  observation  (Fig.  5.6)  shows  a  weak  shear  zone  in  the 
lower  kilometer  of  the  atmosphere  with  the  peak  wind  of  3  m/s  at  200  m  drop- 
ping off  to  less  than  0.5  m/s  at  800  m.   Although  these  conditions  are  not 
particularly  strong,  the  basic  conditions  for  shear-generated  instability  are 
concluded  to  be  present.   The  distance  from  the  upper-air  observation  and  the 
tower  might  also  reduce  the  representativeness  of  the  rawinsonde  observations. 

->.   2        2 
Figure  5.7  shows  the  time  series  of  Ri  ,  9p/9t,  (9U/9z)  ,  and  N   .   The 

time  series  of  9p/9t  documents  the  wave  activity  well,  and  the  acoustic 
sounder  records  in  Figs.  5.8a,b  provide  another  look  at  the  height  distri- 
bution of  the  wave  activity. 

Yet  the  only  clear  evidence  of  wave  propagation  in  the  Richardson  number 

begins  at  2220  and  ends  at  2330  MST.   It  is  interesting  to  note  that  whereas 

■+■    2  2 

the  shear  term  (9U/9z)   appears  to  be  modulated  by  the  wave,  the  N   term  does 

not.   It  appears  that  the  dominant  factor  in  Richardson  number  reduction  is 

the  wind  shear.   Further  investigation  by  decomposing  the  shear  into  its 

azimuthal  and  speed  parts  (Fig.  5.9)  and  comparing  with  Ri,  shows  that  early 

in  the  evening  the  azimuthal  shear  is  large  and  seeks  to  drive  the  Richardson 

number  lower  (Fig.  5.10).   Later,  after  the  Richardson  number  falls  below  0.25, 

the  speed  shear  increases. 

If  we  assume  that  the  onset  of  turbulence  begins  at  Ri  <  0.25,  then  the 
increase  in  wind  speed  might  be  a  result  of  the  momentum  transfer  into  the 
lower  boundary  layer  by  the  turbulence,  while  the  azimuthal  shear  reflects  the 
ageostrophic  conditions  that  created  the  imbalance. 
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Figure  5.5. — Rawinsonde  observations  taken  by  the  National  Weather 
Service  at  Stapleton  Airport,  Denver,  Colo. ,  at  (a)  1700  MST,  6  Novem- 
ber 1981,  and  (b)  0500  MST,  7  November  1981. 
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Figure  5.6. — Wind  speed  vs.  height  for  rawinsonde  observations  at 
1700  MST,  6  November  1981  (solid  line),  and  0500  MST,  7  November  1981 
(dashed  line). 


The  lack  of  modulation  in  N   might  be  an  instrumentation  problem. 

During  calibrations  the  slow-response  temperature  probes  were  found  to  respond 

at  differing  rates  to  changes  in  ambient  air  temperature.   This  will  result  in 

smoothing  the  potential  temperature  gradient.   The  poor  correlation  between 

3p/3t  and  Ri   could  be  due  to  this  factor.   But  the  overall  behavior  of  Ri 
v  v 

through  the  period  shows  evidence  that  it  can  be  tied  to  the  inertial  oscil- 
lations of  the  wind  (Dutton,  1976)  on  a  much  longer  period  than  that  of  acous- 
tic-gravity waves  (Fig.  5.10). 
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Figure  5.7. — Time  series  of  Ri  ,  N   ,  and  (8U/8z)   at  75  m,  and  8p/8t 

measured  at  the  ground. 
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Figure  5.9. — A  azimuth  and  A  speed  over  the  50-  to  100-m  level  of  the 
tower. 
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Figure  5.10. — Time  series  of  Ri  on  an  expanded  scale  showing  long- 
period  oscillation  about  the  critical  Richardson  number. 


5.6   CONCLUSIONS 

The  ability  to  calculate  the  gradient  Richardson  number  using  data  from  a 

boundary  layer  tower  is  demonstrated.   It  is  concluded  that  the  shear  term 

plays  an  important  role  in  the  magnitude  of  the  Richardson  number  and  is  in 

fact  modulated  by  wave  activity.   It  is  also  shown  that  the  behavior  of  the 

shear  term  and  the  magnitude  of  Ri  could  be  a  result  of  inertial  oscillations 

v 

in  the  wind  field  in  addition  to  higher  frequency  wave  activity. 
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